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EDITORS’ PAGE 

Happy New Year! It has been our pleasure to reach you through this journal in the 
year 2003 and we look forward to working with and for you in the ensuing year 2004. 
We thank you for encouraging us through your appreciative words for whatever little 
improvement we could achieve in respect to this journal. However, we most humbly 
admit that a lot more remains to be done. 

Education and research in India has many interesting facets. Both these sectors are 
under tremendous pressure. One curious (and possibly controversial) facet is the 
following : Whatever you do to bring about improvement in these sectors (Education 
and Research), the standards go down. We must remind you that this statement is only a 
partial truth and must be considered with due caution. Research activity of an 
individual or group should be monitored rigorously. A measure of the research activity 
comes from the various publications emanating out of their effort. “Effort” could be a 
genuine concern or the result of a feeling of “competitiveness” amongst fellow 
scientists. In the case of latter, all concerned start “publishing to avoid perishing’’ with 
utter disregard for the quality. Does it not call for caution? Please think! Another 
method adopted for improvement in research and in making it more relevant socially is 
“to increase funding in focussed areas If you do this, many of us have seen in the 
past that the “resourceful and/or ambitious” scientists quickly “de-tune and re-tune” 
their research proposals to snatch the funding and go on an equipment buying spree. If 
you cut funds, the “haves” continue to breathe normally as ever but the earlier “have- 
nots” are made to starve more. Such depressing talks often take place amongst 
scientists as if they themselves are in no way responsible for this state of affairs. Does it 
not call for introspection? 

Assured promotional avenues in research organizations and educational institutions 
to ease tension of scientists have led to disastrous results. Who is to be blamed? Who is 
to be answerable? We leave you to draw your own conclusions since our remarks may 
hurt many. 

Now we come to our journal. The Council of the Academy and we in the Editorial 
Board thought that inclusion of “Lead Articles” by newly elected Fellows/Senior 
Scientists would increase the visibility (and citation index) and usefulness of the 
journal. Yes, it did so. Prof. U.R. Rao in his Inaugural Address delivered at the Annual 
Conference of the Academy held at Ahmedabad in October 2003 also recommended 
publishing of authoritative reviews in Academy’s journal (this address is published in 
this issue of the journal). We are gratified to have received his support and affirmation 
of journal’s policy. However, we ourselves are not satisfied. Only a few newly elected 
Fellows/Senior Fellows are contributing. We request their indulgence to give credence 
to the effort of the Academy to improve this journal. 

Once again, A VERY HAPPY NEW YEAR! ! 

Girjesh Govil 
Jai Pal Mittal 
Suresh Chandra 


The views expressed here are solely those of one of the Editors and do not necessarily reflect those of the Academy 
or the Institute where he works. 




LEAD ARTICLES/OVERVIEWS OF NEW DEVELOPMENTS 


Homology modelling in protein structure prediction : 
Epidermal growth factor receptor kinase domain 


Sunil K. Panigrahi AND Gautam R. Desiraju* 


School of Chemistry, University of Hyderabad, Hyderabad 500 046, India 

*To whom correspondence should be addressed. Tel: 040 23010510 (extn. 4828) Fax: 040 23010567. 

e-mail: desiraiu@.uohyd.ernet.in ~ . , 

Received and Accepted October 20, 2003 


Abstract 

The importance of homology modelling in 
predicting the three dimensional structures of 
proteins has been discussed. Efforts in improving 
the predictive quality are outlined. Successful 
implementation of homology modelling has been 
demonstrated in predicting the structure of the 
kinase domain of Epidermal- Growth Factor 
Receptor (EGFR). A comparison between the 
modelled and experimental X-ray crystal 
structure (C a RMSD 1.96 A) shows the accuracy 
of the modelling methods used. 

Keywords : Biomolecular structure/ computer 
simulation/ cancer/ drug design/ crystallography 


Introduction 

Proteins are among the most important 
biological macromolecules and they are 
extremely versatile in their functions. These 
functions are generally determined by the 
three dimensional structures. Till date, the 
experimental structure of perhaps only 1% 
of all proteins whose sequence is known 
has been determined. This is due to 
practical difficulties inherent to the 
experimental methods (X-ray diffraction, 
synchrotron X-ray diffraction, NMR, 
electron diffraction, electron microscopy, 
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fibre diffraction, neutron diffraction, 
fluorescence transfer) that are used to 
elucidate the structure. Predicting the 

structure of an unknown protein to a near 
enough accuracy through modelling is an 
alternative approach to structure 

determination. Protein structure prediction 
relies on a knowledge of experimentally 

determined structures and is an 

amalgamation of a number of strategies. 
(Fig. 1). The effectiveness of such 
predictive methods depends on the recent 
rapid increase in the size of the Protein 



Fig. 1- Protein structure prediction as an amalgamation 
of various disciplines. 

Data Bank (PBD). 1 Fortunately, there has 
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also been a concomitant decrease in the 
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Fig. 2- Total available structures (green) and annually 
deposited structures (red) of proteins in the PDB. 

number of new folds deposited in the PDB, 
making structure prediction easier. (Fig. 3). 
Lastly, improved hardware and software 
have greatly facilitated this exercise. 


Fig. 3~ Total number of unique folds (red) available in 
the PDB with annually deposited new folds 
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(blue). 

There is now a novel platform to assess 
the strength of modelling methods for 
protein structure prediction. The Critical 
Assessment of Structure Prediction 
(CASP) 2 and the Critical Assessment of 
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Fully Automated Structure Prediction 
(CAFASP) 3 exercises are carried out to 
assess the capability of modern predictive 
methods. CASP is organised every two 
years and evaluates and compares predicted 
structures with crystallographically 
determined ones. The results are discussed 
in an online forum called FORCASP. 
CASP experiments were run in 1994, 1996, 
1998, and 2000, and 2002. Similarly 
CAFASP runs a parallel program solely for 
automated servers and metaservers using 
same targets as those in CASP. All 
CAFASP predictions are part of the CASP 
exercise. These results every two years 
provide an estimate of progress by 
attempting to measure in a quantitative way 
the success of many groups on a predefined 
set of structures. Comparative modelling 
results at CASP 5 have been recently 
summarised. 4 The results show that protein 
structure prediction is becoming more 
reliable, but that there still lie greater 
challenges. 5 A number of methods are 
available (Fig. 4). Among these, 
comparative modelling or homology 
modelling has achieved success. 



Fig. 4- Classification of computational approaches in 
protein structure prediction. 


The present article attempts to explain 
the technical details of this method and 
illustrates these with respect to the 
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modelling of the Epidermal growth factor 
receptor (EGFR) kinase domain. EGFR has 
attracted attention because of the findings 
that deregulation of this receptor system is 
a significant factor in the genesis or 
progression of several human cancers. 6 
Shortly after we completed our modelling 
study, the experimental crystal structure of 
EGFR was published. 7 This gave us an 
opportunity to compare our model with the 
experimental structure. 

A review of homology modelling 

In homology (comparative) modelling 
the three dimensional structure of a protein 
(target) is built based on its alignment to 
proteins of known structure (templates). 
The method has been reviewed extensively 
in the literature. 8 The basic assumptions 
and steps are as follows: 

• Related members of a family have 
evolved divergently with respect to 
sequence and structure. 

• Two evolutionary related neighbours 
with a similar primary sequence will 
have similar three dimensional 
structures. 

• Two sequences are evolutionary related 
if the similarity between the sequences 
is greater than 30%. The model is best 
when they are >50% identical, and not 
very useful when <30%. 

• The target sequence is used in a 
sequence similarity search of the PDB 
to determine if it has recognizable 
similarity to sequences for which a 
structure is known. 

• Sequence alignments are then 
constructed between the sequence to be 
modelled and the template sequences. 
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• Multiple domain containing proteins are 
modelled by considering individual 
domains separately. 

Steps involved in homology modelling 

The various steps involved in homology 
modelling are illustrated in Fig. 5. These 
are now detailed. 



Fig. 5-Flow chart showing the steps involved in 
homology modelling. 

1. Selection of the template 

This exercise uses various search 
methods and database mining tools to 
identify templates in databases such as 
PDB, 1 SCOP, 9 DALI, 10 and CATH. 11 
Broadly protein comparison methods can 
be classified into three main classes. 

a. Pairwise sequence-sequence compa- 
rison. These methods compare the 
target sequence with each of the 

3 
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database sequences independently, 
using pairwise sequence-sequence 
comparison (FASTA 12 and BLAST 13 ). 

b. Multiple sequence comparison methods. 
These rely on multiple sequence 
comparisons to improve the sensitivity 
of the search. 14 For a given sequence, 
an initial set of homologues from a 
sequence database is collected, a 
weighted multiple alignment is made 
from the query sequence and its 
homologues, a position-specific scoring 
matrix is constructed from the 
alignment, and the matrix is used to 
search the database for new 
homologues. These steps are repeated 
until no new homologues are found 
(PSI-BLAST 14 ). 

c. Threading or 3D template matching. 15 
These methods rely on pairwise 
comparison of a protein sequence and a 
protein of known structure. The target 
sequence is threaded through a library 
of 3D folds. For each sequence- 
structure pair, independently, structure- 
dependent scoring functions are used to 
judge whether or not a given target 
sequence adopts any one of the many 
known 3D folds as predicted by an 
optimization of the alignment. 

2. Alignment of target and template 

The next step in comparative modelling 
is the alignment between the target 
sequence and the template structure. 
Alignment means a rearrangement of 
subsets (of a sequence) to find maximum 
similarities between the target sequence and 
the template structure. 16 It consists of 
establishing equivalences between the 
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residues in two different proteins, as is the 
case with conventional sequence alignment. 
The alignment is done on the basis of the 
three-dimensional coordinates correspond- 
ing to each residue, and not on the basis of 
the amino acid type of the residue. 

The procedure involves the alignment 
of potential templates, which are prepared 
by superposing their structures. Next, the 
sequences that are clearly related to the 
templates and are easily aligned with them 
are added to the alignment. The same is 
done for the target sequence. Finally, the 
two profiles are aligned with each other, 
taking structural information into account 
as much as possible. 17 In principle, most 
sequence alignment and structure com- 
parison methods can be used for these 
tasks. Widely used is the alignment based 
on dynamic programming. Their success 
depends on the degree of sequence identity 
between the target and template. For 
closely related protein sequences with 
identity over 40%, alignment is almost 
always correct. Regions of low local 
sequence similarity become common when 
the overall sequence identity is under 
40%. 18 Alignment becomes difficult in the 
“twilight zone” of less than 30% sequence 
identity. 19 As the sequence similarity 
decreases, alignments contain an increa- 
singly large number of gaps and alignment 
errors, regardless of whether they are 
prepared automatically or manually. 
Because insertion and deletion occur at the 
loop regions of the protein between 
secondary structures, improvements can be 
made by introducing penalties for 
insertions/deletions of a-helices or (3- 
sheets. Phylogenetic inference and motif 
comparison in the alignment step are vital 
for correct modelling. 

SCIENCE LETTERS, VOL. 27, NO. 1 & 2, 2004 
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3. Model building 

After alignment the model is built using 
various model building methods. The test 
of any modelling method lies in its 
accuracy when the sequence similarity is 
less than 40%. Most important of all is the 
selection of the template and proper 
alignment, which is crucial for model 
quality irrespective of the modelling 
methods chosen. Three model building 
methods are generally used: 

(a) Model building by rigid body assembly 
method. 20 The method dissects the 
overall protein folds into structurally 
conserved regions, variable loops, and 
side chains followed by model building 
with a small number of rigid bodies 
obtained from aligned protein 
structures. 

(b) Segment matching method, 21, 22 relies 
on the approximate positions of 
conserved atoms in the templates. 
Comparative models are constructed by 
using a subset of atomic positions from 
template structures as “guiding” 
positions, and by identifying and 
assembling short, all-atom segments 
that fit these guiding positions. The 
segment matching method is based on 
the idea that most hexapeptide segments 
in proteins can be clustered into 
approximately 100 structural classes. 23 
Therefore, by transferring the 
coordinates from the template segment- 
wise; it is possible to generate a model 
which has a segment structure very 
close to the actual one. 

(c) Model building by spatial restrain 
satisfaction, 24,25 uses either distance 
geometry or optimization techniques to 
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satisfy spatial restraints obtained from 
the alignment. These methods generate 
many constraints or restraints on the 
structure of the target sequence, based 
on its alignment to related protein 
structures. The restraints are generally 
obtained by assuming that the 
corresponding distances and angles 
between aligned residues in the 
template and the target structures are 
similar. These homology derived 
restraints are usually supplemented by 
stereochemical restraints on bond 
lengths, bond angles, dihedral angles, 
and nonbonded atom-atom contacts 
obtained from a molecular mechanics 
force field. The model is then derived 
by minimizing the violations of all the 
restraints. 

In the process of model building the 
loop as well as the side chain is built. 
Further improvement in the quality of the 
loop and the side chains may be obtained 
through specialised modelling tools. 

4. Loop building 

Loop modelling methods can be 
categorized into ab initio , database search 
and combinations of both methods. 

a. Ab initio methods : 

De novo design of the loops is carried 
out based on the conformational analysis 
guided by a scoring or energy function in a 
particular environment. 46, 47 Many methods 
exploit different protein representations, 
energy function terms, and optimization or 
enumeration algorithms. All methods 
follow basic steps of conformation 
generation, sampling of the conformation 

5 
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and selection of the appropriate confor- 
mation. 

b. Database search techniques 26 ' 29 

Loop prediction by a database approach 
finds a segment of the main chain that fits 
the two stem regions of a loop. The stems 
are defined as the main chain atoms that 
precede and follow the loop but are not part 
of it. The database approaches are best 
suited for proteins where loops have a 
small conformational space. The 
exponential increase in the number of 
geometrically possible conformations as a 
function of loop length limits the database 
approach. Another approach of loop 
modelling follows a combination of both ab 
initio and database search methods 
implemented simultaneously. 

5. Side chain modelling 

The side chains play a major role in the 
overall packing quality, interactions types 
and stability of the protein. Side chain 
modeling is a daunting task. Some methods 


use energy functions to predict side chain 
conformations. 30 Prediction of side chain 
conformations are done along with 
backbone prediction. Side chain prediction 
at the structurally conserved region is easier 
than in the structurally variable regions. 
Inclusion of solvation and hydrogen 
bonding terms improves the success of side 
chain prediction markedly. 31 There is a 
rapid decrease in the conservation of side 
chain packing as the sequence identity falls 
below 30%. 32 Although the fold is 
maintained, the pattern of side chain 
interactions is generally lost in this twilight 
range. Various approaches to side chain 
modelling have been reviewed elsewhere. 33 

6. Model evaluation 

Quality control is foremost in model 
building. The built model should satisfy (a) 
Stereochemical accuracy, (b) Packing 
quality and (c) Folding reliability. Various 
quality checks are listed in Fig. 5. Protein 
quality check tools like PROCHECK, 34 
WHATCHECK, 35 AQUA, 36 SQUID 37 are 
available to judge . protein quality. 



Fig. 6- Protein quality check parameters. 
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7. Model Refinement 

Refinement is done by further 
minimisation to eliminate bad contacts. 
Sometimes the model is subjected to 
molecular dynamics simulations to find out 
the lowest energy conformation. A manual 
adjustment of the packing quality of side 
chain improves model quality. 

Errors in homology modelling 

A model quality falls with the decrease 
in the similarity between the target and the 
templates. Frequently a model may go 
wrong if (a) Incorrect templates are 
selected, especially when distantly related 
proteins are used as templates (i.e., less 
than 25% sequence identity), (b) The 
alignment is not proper. If the target- 
template sequence identity decreases below 
30% then the model accuracy depend solely 
on alignment quality. Alignment errors can 
be reduced by constructing a multiple 
alignment taking a large number of 
sequences. Multiple alignments are 
generally more reliable than pairwise 
alignments. However in cases where the 
sequence identity is greater, the pairwise 
alignment also performs well. 


Strengths and lacunae of homology 
modelling 

A model may attain a structure close to 
the experimental one when the template 
selection is good (family, similarity, 
tertiary structure). Then, methods to 
construct conserved and variable regions 
and the side chain work properly. Model 
prediction fails when the target and 
templates have a large number of insertions 
and deletions so that the alignment results 


in improper coordinate transfer, when the 
target and template have sequence 
homology less than 30% resulting in 
improper alignment, and when templates 
having poor resolutions are chosen. When 
the prediction is not so good, the final 
model may possess more of template 
character than target character. To 
summarise, the different methods in 
predicting three dimensional structure are 
illustrated in Fig. 7. 



Fig. 7-Overview of the method as a funotion of 
sequence similarity to template (expressed as a 
percentage). 


Results and Discussion 


EGFR belongs to the larger class of 
trans-membrane growth factor receptor 
tyrosine kinases (RTKs). 38 It consists of 
1 186 amino acid residues out of which the 
catalytic domain contains 310 residues. 39 
The receptor is divided into three main 
regions: an extracellular domain ligand 
binding part, a hydrophobic membrane 
spanning region, and a cytoplasmic portion 
containing the kinase domain.(Fig. 8). 
EGFR and its ligands (EGF, TGF-a) have 
been implicated in numerous tumours of 
epithelial origin (e.g. squamous cell 
carcinoma; breast, ovarian, NSC lung 
cancer etc) and proliferative disorders of 
the epidermis such as psoriasis. 6 
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A number of competitive inhibitors 
targeting the ATP binding site are in 
various phases of clinical trials. In an 
ongoing research in our group in 
Hyderabad we wanted to study the nature 
of the interaction between 4-anilinoquin- 
azolines, an important class of inhibitor, 
and the kinase domain. 40 Because the 
crystal structure was not available at the 
time when we carried out this study, we 
built the homology model of the kinase 
domain. Modelling was carried out using a 
Silicon Graphics 02 workstation and the 
‘Homology’ module of MOE (Molecular 
Operating Environment, Chemical Compu- 
tation Inc). 41 While our study was nearing 
completion, the crystal structure of the 
protein was published. 7 Thus a comparison 
of structure quality was performed between 
our model and the experimental crystal 
structure. 


Methodology 

1. Template selectipon 

A total of 263 amino acids (688 to 950 
of the original receptor) out of 3 1 0 residues 
of the kinase domain ware taken from 
SWTSS-PROT, accession no P00533. 
Modelling was carried out selecting 
IR3A.pdb 42 as the primary template. Both 
share a sequence similarity of 35%. To 
enhance the quality of proteins in terms of 
the structurally conserved regions and 
structurally variable regions we selected 
3LCK.pdb, 43 lIRK.pdb, 44 lFPU.pdb, 45 
lCDK.pdb 46 as secondary templates. The 
criteria of selecting IR3 as the primary 
template are as follows: 

a. It is a membrane receptor having 
architecture similar to that of EGFR. 
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b. The signal transduction pathway for 
insulin receptor is quite similar to 
EGFR . 

c. The activated form of the insulin 
receptor has the same type of 
conformation in the catalytic domain of 
the activated form of EGFR. 

2. Alignment 

The target and the templates were 
aligned using MOE-Align module which is 
a modified version of the Needleman and 
Wunsch 47 method. Alignment was compu- 
ted by optimising the function based on 
residue similarity score obtained from the 
Gonnet 48 matrix method using a gap of one 
and gap extension of one. Four automated 
steps, namely initial pairwise build-up 
(Progressive/Tree-based), round-robin rea- 
lignment, randomised iterative refinement, 
and structure-based realignment of MOE- 
Align module were used for the alignment 
of target and templates. 

3. Model building 

Model building was carried using the . 
following steps: 

(a) Initial partial geometry specification. 
Initial geometry from regions of all 
template chains were copied. All 
coordinates were copied where residue 
identity was conserved. 

(b) Intermediate model construction. Using 
a Boltzmann-weighted randomized 
modelling procedure combined with 
specialized logic for proper handling of 
insertion and deletions, 49 independent 
models of the target proteins were built. 
Models were evaluated for residue 
packing quality. Automated loop and 
side chains modelling was selected. 

SCIENCE LETTERS, VOL. 27, NO. 1 & 2, 2004 




Fig. 9- Superposition of our modle (original colours) with the experimental crystal 
structure (blue) of EGFR kinase domain. 
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(c) The final model was built and 
minimized. Coordinates were generated 
as the average of the atom coordinates 
of the intermediate models. After the 
final model was built, a rotamer 
exploration was carried out for 
important residues in the catalytic 
domain. The lowest energy 
conformation was chosen for the final 
model. 

4. Quality check (evaluation, refinement and 
validation) 

The final model was subjected to 
quality check. MOE stereochemical quality 
evaluation tools were applied to confirm 
that the stereochemistry of the model is 
reasonably consistent with typical values 
found in crystal structures. Tools like 
PROCHECK and Ramachandran plot were 
used to evaluate the model quality from a 
stereochemical viewpoint. Further energy 
minimisation was carried out to reduce bad 
contacts. The catalytic domain was treated 
with special care. Docking studies were 
performed with one of the potent inhibitors 
to validate the quality of the model. 

Comparison with the crystal structure 

The superposition of the homology built 
model and the experimental crystal 
structure M14.pdb is shown in Fig. 9. The 
RMSD with respect to the C a atoms was 
1.96 A. Both show the bilobate nature like 
other kinases. The amino terminal lobe 
consists mostly of p -strands and one a- 
helix (aC) whereas the carboxy terminal 
lobe consists mostly of a-helices. The 
binding site of ATP, ATP analogues and 
ATP -competitive inhibitor is a cleft created 
by the amino and carboxy terminal lobe. 
This cleft consists of five important 
regions, the glycine-rich nucleotide 
phosphate-binding loop (Gly695-Gly700) 
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of the amino terminal lobe, the DGF motif 
(Asp831-Gly833), the presumptive 
catalytic (general base) Asp813, the 
catalytic loop (Arg812-Asn818) and the A- 
loop (Asp831-Val852) of the carboxy 
terminal lobe. These regions are indicated 
in Fig. 9. A comparative look at different 
regions of the model and the crystal 
structure gave us greater confidence on the 
accuracy of the model quality. 

Conclusions 

Homology modelling is a successful 
means of protein structure prediction. We 
have successfully implemented this method 
to build the three dimensional structure of 
epidermal growth factor receptor kinase 
domain. The model quality is as good as the 
reported crystal structure of the protein. 
However homology modelling is not the 
only method to predict the structure of a 
protein. Threading and ab initio methods 
have shown promising results in model 
building when the percentage, similarity is 
less than 30% between template and target. 
In threading or fold recognition, a library of 
unique or representative structures is 
searched for structure analogues to the 
target sequence, and is based on the theory 
that there may be only a limited number of 
distinct protein folds. 50 Modelling of the 
protein when there is no structural template 
is existence is quite difficult. However ab 
initio approaches or de novo design based 
on biophysical principles have gained some 
success in this regard. Ab initio methods 
attempt to predict the structure from 
sequence alone, without relying on 
similarity at the fold level between the 
modeled sequence and any known 
structure. 51 It attempts to find the global 
free energy minimum by an exploration of 
many conceivable protein conformations. 
These methods assume that the native 
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structure corresponds to the global free 
energy minimum accessible. 

All in all, modelling methods have 
achieved some success in predicting the 
structures of unknown proteins but there 
still remain many challenges to overcome. 
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Abstract 

An overview of various bivariate geometric 
distributions is given. The focus is on 
characterization results. The role of bivariate 
distributions in reliability and Statistical Quality 
Control is illustrated, with examples, in detail. 

(Keywords : characterization/survival function/ 
failure rate function/mean residual life function/ 
lack of memory property). 

Introduction 

The extension of univariate probability 
distributions to two or more dimensions is 
not an easy task. Very often investigations 
in this regard result in more than one 
distribution. Researchers generally believe 
that certain basic properties associated with 
a univariate distribution must be preserved 
(in an appropriate form) by their counter 
parts in higher dimensions. One of the main 
problems confronted in pursuing this 
approach is the fact that there may not be a 
unique formulation of such a basic property 
of a univariate distribution in higher 
dimensions. Also, researchers consider in 
higher dimension, real life situations 
similar to those in one dimension for which 
a univariate distribution is known to be a 
suitable stochastic model and investigate 
stochastic models appropriate for those 
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situations. Naturally there will be several 
candidate probability distributions in higher 
dimensions which could be considered as 
extensions of a univariate probability 
distribution. 

The literature concerning bivariate and 
multivariate exponential distributions is 
very rich (see, for example, Galambos and 
Kotz 1 ). The investigations were motivated 
by consideration of appropriate stochastic 
models as well as by consideration of 
appropriate versions of characterizing 
properties of univariate exponential 
distribution. However, comparatively less 
is known about multivariate version of the 
geometric distribution, even though it plays 
an important role as the discrete counterpart 
of the exponential distribution in stochastic 
modelling. It will be seen that different 
ways of looking at a problem resulted in 
different versions of geometric distributions 
in higher dimensions, even in the modest 
literature concerning this. Interestingly, 
some researchers have formulated a 
bivariate geometric distribution whose 
marginals are not geometric. 

Bivariate geometric model was initially 
discussed by Downton 2 in connection with 
consideration of shock models in reliability 
that lead to bivariate exponential 
distribution. Commenting on Downton’ s 
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model, Hawkes 3 formulated a bivariate 
geometric distribution through bivariate 
Bernoulli distribution. Around the same 
time, Paulson and Uppuluri 4 , Esary and 
Marshall 5 and Arnold 6 discussed the 
bivariate and multivariate geometric 
distributions from the point of view of 
compounding models. 

Block et a/. 7 gave a very general 
definition of a bivariate geometric 
distribution as one for which the marginals 
are geometric and gave several examples. 
They constructed some auto regressive 
processes with bivariate geometric 
marginals. 

Another approach adopted in the 
literature concerning bivariate geometric 
distribution is by consideration of 
characterizing properties. 


R{m,n) = 

f p" p”~ n 0 <n<m 
[ P m Pi "" 0 <m <n , m,n = 0,1,2,.. 

where 0 <p <pt <1, /— 1 ,2 ,p\+p 2 < 1 +p, 

or equivalently by its probability, mass 
function (pmf) 

p" pr~ x ( A - P)( l - Pi ) if 0 < « < m 
• p m (l + p~ Pi ~ p 2 ) if 0 < m = « 

P m P2~ m ~ ] (Pi ~ p)Q~ p 2 )tf0< m <n 

m,n = 0,1,2,... 

where 0 <p<p» i- 1,2 and p x + p 2 < 1 +p. 


The aim of this paper is to review the 
literature concerning characterization 
results for several bivariate geometric 
distributions. We shall also mention 
relevant results concerning a multivariate 
geometric distribution. No proofs are given 
but adequate references are provided. 


Remark : If p x + p 2 = ^+P then the 
distribution has no mass at the points 
m=0,l,2„... 

(b) The BGD introduced by Phatak and 
Sreehari 8 , henceforth referred to as BGD(P 
& S), is given by the pmf 


The outline of this paper is as follows. 
We first list various bivariate geometric 
distributions. We then discuss some 
definitions and concepts that play important 
role in characterizing univariate geometric 
distribution. Finally we discuss, without 
proofs, the characterization results 
concerning the various bivariate geometric 
distributions listed earlier. 

Various bivariate geometric distributions. 

(a) The bivariate geometric distribution 
introduced by Hawkes 3 , henceforth referred 
to as BGD(H) is given by its joint survival 
function 


f{m,n) = 




v « j 


Pi P" 2 (1 ~ Pi ~ Pi ),m,n = 0,1,2,. 


where 0 <p x ,p 2 < p x + p 2 < 1. 

(c) Nair and Nair 9 introduced a BGD, 
henceforth referred to as BGD(N & N), by 
its joint survival function 

R(m,n) = p x pi 9 mn ,m,n = 0,1,2,... 
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where 

0 < p x ,p 2 < 1,0 < 0 < 1,1 — 0 < (l-p,0)(l-/? 2 0) 

Corresponding pmf is given by 

fim,n) = P : P n 2 e mn - x [(\-p,Q n+ ')(l-p 2 0 m+I ) 

+ 6 - 1 ] 

where p x , p 2 and 0 satisfy the above 
inequalities. 

(d) Nair and Asha 10 considered a BGD, 
henceforth referred to as BGD(N&A), 
which is given by the joint survival 
function 

R(m,ri)= 

f p: if p:^p2 

IP 2 f Pi<P\ » rn,n = 0,1,2,- 

where 0 < p, <1, z- 1,2.... 

(e) Asha and Sankaran 11 introduced a 
BGD, henceforth referred to as BGD(A & 
S), which is given by its pmf 

ffii,n)= 

(Pi^ 2 )" ( P 2 P 2 0 - P\ )0 ~ p \ ) f 0 ^ m < n 

' (PiPif ( 1 -Pi) 0 -/ J 2 ) i/0<m = n 

(P\P2 )" Oi* r'"' 1 Pi (1 - Pi ) 0 - P\ ) if 0 Z n < rn 

m,n = 0,1,2,... 

where 0 < p t , p] < 1, p ’ * p x p 2 , z= 1 ,2. 

Remark : All the bivariate geometric 
distributions except the last have geometric 
marginal distributions. 


Univariate geometric distribution 

We shall now state some definitions and 
discuss related basic properties of a 
univariate geometric distribution with 
support N = {0,1,2...}. We mention a real 
life situation for which a geometric 
distribution is an appropriate stochastic 
model. Unless stated otherwise all random 
variables (r.v.s) are nondegenerate 
nonnegative integer valued. 

Definition 1 : A r.v. X is said to have a 
geometric distribution if 

P(X=n)= p(\- p y ( 1 ) 

n — 0,1,2,..., 0 <p <1. 

We shall then write X~G(p). 

Definition 2 : The failure rate function 
of a r.v. X is given by 

h(n) = P(X= n\X>n)= P(X= n) / R(n) 

where R(n) is the survival function, P(X > 
n). 

Whenever P(A\B ) or E{Z\B) is 
considered for events A,B and r.v. Z, it is 
assumed that P(B) > 0. 

Definition 3 : Residual life of a r.v. X 
after an elapsed time x is denoted by 

{X-x \X>,x} and the mean residual life 
function p(x) is given by 

ju(x) = E[X-x\X>x]. 

Definition 4 : A r.v. X is said to possess 
the lack of memory property (LMP) at t if 
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P[X> t+s \X> t] = P(X> s), s > 0. 

We need the following results 
concerning a. geometric r.v. for later 
developments and they are easily checked. 


If X~G(p) then 

R{n) = (\-p) n ,n = 0,1,2,... (2) 

h (ri) = l-p, n = 0,1,2,... (3) 

|u(«) = l+£(AD,n = 0,l,2... (4) 

and 

if has LMP at alU = 1,2,. . . (5) 


It is well known that each of the results 
(2) - (5) uniquely determines the geometric 
distribution in the class of discrete 
distributions. In fact, if X has LMP at t =1 , 
then also X has a geometric distribution. It 
is of interest to note that X need not have 
geometric distribution if if has LMP at only 
one point t > 1 . (see Phatak and Sreehari 8 , 
Sreehari 12 and Srivastava and Bagchi 13 ). 
However, if X has LMP at two points m and 
n which are relatively prime then X has 
geometric distribution (Sreehari and 
Srivastava and Bagchi 13 ). In fact, 
Srivastava and Bagchi proved a slightly 
more general result. 

We shall now discuss some properties 
of a geometric distributuion whose 
analogues in higher dimensions will be 
discussed later in the paper. 

Property 1 

In the class of power series distributions 
given by P (X = n) cc a„Q ", n =0,1,2,..., 

16 


geometric is the only one with E (Jf)=0 /(l- 

0). 

Property 2 

If X and Y are independent identically 
distributed r.v.s then 

P (X+Y= k) oc (k+l) P (X= £), £=0,1,2,... (6) 

if and only if their common distribution 
is geometric. 

Property 3 

For some r.v. Y 

P (Y < n) - P (X+ Y< n) oc P (X+ Y = n), 
n = 0,1,2,... (7) 

if and only if X is a geometric r.v. 

If X~G(p), Y=X+l takes values 1,2,... 
and is called a modified or truncated (at 
zero) geometric r.v. 

Example. 

Consider a system which is subjected to 
a series of shocks each of which has a 
probability p of making the system 
‘nonworking’. Suppose that the shocks 
occur independent of each other and that no 
damage to the life of the system occurs by a 
shock which is not ‘fatal’. Then the number 
X of non-fatal shocks received by the 
system before it received a shock that is 
‘fatal’, has the geometric distribution given 

at (1). 
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Bivariate geometric distribution, 
BGD(H). 

In this section we shall discuss an 
example which leads to BGD(H) and give a 
number of characterization results 
concerning this distribution. In this and 
later sections all bivariate random vectors 
have support N x N where N = {0,1,2,...}, 
unless otherwise stated. 

Example : 

Suppose two components Ci and C 2 of 
a system receive shocks which occur in 
independent Poisson streams of rates A,i and 
X 2 . Suppose further that each shock* may 
not necessarily be fatal resulting in the 
failure of the component. Let X, be the 
number of shocks, received by the 
component C, before its failure. Suppose 
the joint probability generating function 
G(si,s 2 ) of (X U X 2 ) is 

G(s\, s 2 ) =sis 2 / (1 +a+p+y - a - 0 s 2 - y 
s\ s 2 ) 

where 0< a, 0, y<l. 

We note that the marginal distributions 
of X\ and X 2 are geometric distribution over 
{1,2,...}. Using the bivariate geometric 
distribution given by Downton 2 

discussed the joint distribution of the life 
times of the two components Q and C 2 
involving the five parameters Xi, X 2 , a, 0 
and y. Commenting on Downton’ s work, 
Hawkes 3 introduced a bivariate geometric 
distribution, through a bivariate Bernoulli 
distribution, with support {1,2,3,...} x 
{1,2,3,...}. This distribution is one of the 
widely discussed bivariate geometric 
distributions and is a reformulation of 
BGD(H). 
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The illustration given in the previous 
section can be reformulated in terms of 
observation of a sequence of independent 
Bernoulli r.v.s {X n }, for which P (X„ = 1) = 
p and P ( X n = 0)=l-p, till a ‘one’ is 
observed. Our example of BGD(H) is given 
in terms of a bivariate Bernoulli 
distribution. 

Definition 5 : Vector ( X,Y) is said to 
have a bivariate Bernoulli distribution with 
parameters ( p 00 , p 0l , p lQ ) if 

P(X= 0, Y= 0) =p 00 

P(X= 0, T= 1) = p 0l 

p(x= i, r= o) = Pl0 

P(X— 1, Y— 1) l-p 0 o ~Po\ ~ P\q ~ Pw • 

For convenience, let us denote 

Pi o +Pn =1 P«+ andpo,- +pij=p +j , ij = 0,1. 

Consider a sequence {(X„;Yn)} of inde- 
pendent identically distributed bivariate 
Bernoulli random vectors. Let U denote the 
number of 0’s before the first 1 in the 
sequence X\, X 2 ,..., and let V denote the 
number of 0’s before the first 1 in the 
sequence Y\,Y 2 ,... Then the joint 
distribution of JJ and V is given by 

P(U=m, V=n) = 


„/w 

Poo 

n-m-l 

P + 0 

-PlO P+1 

if m<n 

Poo 

p ’ll 


if m-n 

Poo 

OT-A7-1 

P o + 

Pol Pi + 

if n<m 
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Reparameterizing this distribution we 
note that (U,V) has the BGD(H). To see 
this, take p=p 00 , p x = p 0+ and p 2 =p +Q . 


where 0< Cj < 1, i=T,2,3 and the marginal 
distributions of U and V are geometric with 
parameters 1 -p l and 1 - p 2 respectively. 


Marshall and Olkin 14 observed that — 

4 

< corr ( U,V ) < 1 . 

Further, the marginal distributions of U 
and V are geometric ( 1 - p, ) and geometric 
(1 -p 2 ) respectively. 

A natural extension of the failure rate 
function to two dimensions is the bivariate 
failure rate function 


Remark : It is noted that the proof given 
by Asha and Nair is not consistent with 
their notation for the survival function. 

Another extension of the failure 
(hazard) rate function to two dimensions is 
by considering the vector failure rate 
function ( h\{m,n ) , h 2 (m,n)) where 

h\(m,ri)=P(U=m , V> ri)/R(m,ri) and 


h(m,n ) = 


P(£/ = m,F = w) 
P(U >m,V>n) 


This was also called Scalar failure rate 
function by Asha and Nair 15 . It is easily 
seen that h(m,n ) is not a constant for the 
BGD(H). In fact, it is constant in each of 
the regions 0 <m <n, 0 < m = n and 0 < n 
< m. Asha and Nair proved the following 
result which gives a characterization of the 
BGD(H). 

Theorem 1 : A random vector ( U,V) has 
the BGD with survival function 


R(m,n)= 


P Pi 

_ n-m 

P P 2 


n<m 

m<n 


( 8 ) 


0 <p< Pi < \,j=l,2,p\ +p 2 < 1 +p if and only 
if the Scalar failure rate function is of the 
form 


h 2 (m,n)=P(U > m, V~ n )/ R(m,n), 

whenever R ( m,n ) >0. In terms of this Nair 
and Asha 16 obtained the following. 

Theorem 2 : A random vector ( U,V) has 
the BGD(H) if and only if the components 
of the vector failure rate function are of the 
forms 


h\{m,n) = 


a, 


for n<m 


a 2 for m<n 


h 2 (m,n ) = 


for n<m 
I b 2 for m<n 


for all m,n = 0 , 1 , 2 ,... and some constants 
0 <a i} bj < 1 , i=l ,2 satisfying ( 1 - a\) ( 1 - b\) 
= ( 1 - ai) (1- bp) and 2>i(l- a\) < 62 . Here p\ = 
1 - a\, p 2 = 1 -T >2 and/?=(l- ai) ( 1 - b\). 


h(m,n ) = 




if n<m 
if m = n 
if m<n 


Remark : The above result can be 
reformulated in terms of vector mean 
residual life time, as was done by Nair and 
Asha. 
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If ( U,V) has the BGD(H) it is easily 
verified that the conditional distribution of 
U given V=n is given by its pmf g(m \ n) - 
P(U=m | V=n) where 

g(m | n ) = 

P\~ X (P 1 P\Pi)" [0-A)(A -P)f(}-P 2 )] 

if n < m 

, [(1 + P- P\ -Pi)l(\- Pi)} (p/PiY if m = n 
p m p- (m+]) [p 2 -p] if m<n 

It is well known that if the conditional 
distribution _/(*i U 2 ) of A”i given X 2 =x 2 , and 
the marginal distribution of X 2 are known 
then the joint distribution of (X\, X 2 ) is 
determined. But knowledge of/(*i I * 2 ) and 
that of the marginal distribution of X\ is not 
adequate to determine the joint distribution 
of (Xu X 2 ). In this context Nair and Asha 16 
proved an interesting characterization result 
based on a distributional property. Similar 
results for BGD(P&S) and BGD(N&N) are 
given in Theorems 12 and 15. 

Theorem 3 : If a random vector (U,V) 
has P(U=m \ V=n) given as above then U is 
geometric (l-p\) if and only if V is 
geometric (\-p 2 )- Consequently, (U,V) is 
BGD(H). 

It is well known that if X\, X 2 ,...Xn are 
independent geometric r.v.s then min(Xi, 
X 2 ,...Xn) is a geometric r.v. Similarly, if 
(X u Y x ), (X 2 ,Y 2 )...,(X n ,Y n ) are independent 
random vectors having bivariate geometric 
distributions of the form given by (8) then 

f } 

min X„ min Y, 

U <i<n 1 <i<n J 
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has a BGD(H). 

An interesting result connecting the 
BGD(H) with minima of independent 
geometric r.v.s is the following. 

Let Y\, Y 2 , Yj be independent geometric 
r.v.s with parameters (1 -a), (1-/3) and (1 -y) 
respectively. Then (U,V) where 
L=min(Yi,F 3 ) and V=min(Y 2 ,Y 3 ) has 
BGD(H) given at (8) with p\ = ay, p 2 = py 
and p-afiy. This interesting result noted by 
Block et al. 7 is a part of Theorem 5.1 in 
Nair and Asha 16 and is useful in generating 
random vectors having BGD(H) given at 
(8) for which p\p 2 < p. 

Remark : Nair and Asha 16 claim a 
converse to this result in Theorem 5.1 of 
their paper which appears to be incorrect. 

Nair and Asha considered the bivariate 
no-aging property 

P[X x > m+k, X 2 > n+k | Xi> k, X 2 > k] = P[X x > 
m, X 2 > n] (9) 

for m,n,k= 0,1,2,... to investigate the 
behaviour of the joint distribution of (X x , 
X 2 ). It is noted that this is a kind of bivariate 
LMP. They proved the following: 

Theorem 4. A random vector {U,V) has 
the bivariate no-aging property given at (9) 
and the marginal distributions are 
geometric if and only if the joint 
distribution of (U,V) is BGD with survival 
function R(m,n ) given in Theorem 1 for 
some p x ,p 2 and p. 

Failure of a two-component parallel 
system can be considered to consist of two 
stages: first one of the two component fails 
and later the remaining one fails. Working 
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on this, Sun and Basu 17 introduced total 
failure rate for the discrete r.v.s set up. 

Definition 6 : For a random vector(£/,F) 
the vector ( r(t ), rfm | n), r 2 (n | m)) is called 
the total failure rate where 

r(t) = P [min(£/,F) = t\u>t,V>t] (10) 

rfm | n) = P [U=m, V=n]/P ( U> m, V=n) for 
m>n (11) 

r 2 (n \m) = P [U=m, V=n]/P (U=m, V>n) 
ioxn>m. (12) 

Note that r(f) is the failure rate of 
min (U,V) whereas rfm | n) is the 
conditional failure rate of U given XJ>V and 
V- n. Similar meaning is conveyed by 
r 2 (n | m). 

Sun and Basu obtained the following 
characterization result of BGD(H) in terms 
of total failure rate. 

Theorem 5. A random vector ( U,V) has 
the BGD(H) if and only if it has geometric 
marginals and a constant total failure rate. 
In this case, r(f)= 1 -p, rfm | n) =1 - p x for m 
> n and r 2 {n I m) = 1 - p 2 for n > m. 

Asha and Nair 15 also considered the 
conditional failure rate functions rfm | n) 
and r 2 (n | m) and studied their behaviour in 
the context of BGD(H). 

Remarks : Suppose that a random vector 
(U,V) has the BGD(H) given by (8). Let 
N\=U+1 and N 2 = V+ 1 . The random vector (N\, 
N 2 ) takes values in { 1 , 2 , . . . } x { 1 , 2 , . . . } and its 
distribution may be called modified BGD. This 
distribution has been considered by Downton 2 , 
Paulson and Uppuluri 4 , Esary and Marshall 5 , 
Arnold 6 , Block and Paulson 18 , Marshall and 


Olkin 14,19 , Block et al. 7 , Nair and Nair 20 , Wu 21 , 
Nair and Asha 10 and Abraham and 
Balakrishnan 22 among others. 


Block and Paulson 18 and Block et al. 1 
considered bivariate geometric compound- 


ing, namely 


( N 2 

(5.n)= 2*,. T.r, 

V ' 1 


where 


(X\,Y\), (X 2 J 2 ),... are independent 

identically distributed random vectors, 
(N\ ,N 2 ) has the BGD discussed above and 
(Nu N 2 ) is independent of {{X n , Y„)}. Block 
and Paulson investigated the infinite 
divisibility of the distribution of (£, p) . 


Marshall and Olkin 19 considered the 
stability of geometric maximum and 
geometric minimum. Wu 21 also considered 
bivariate geometric summation and studied 
the limiting distributional properties. Wu 
called the modified BGD of (N\, N 2 ) as the 
general BGD. Wu called the same as 
ordinary BGD in case P(N\ = N 2 )=0, which 
holds if 1 +p ~ P] + p 2 ■ For a discussion on 
a general multivariate geometric 
distribution we refer to Arnold 6 . 


n 

Block et al. constructed bivariate 
exponential (geometric) auto regressive and 
auto regressive moving average models and 
investigated the stationarity of BEAR and 
BGAR. 

Abraham and Balakrishnan considered 
Bivariate exponential auto regressive 
(BEAR(l)) process using the modified 
BGD in the context of a repairable system. 


Bivariate geometric distribution, BGD 
(P & S). 

Phatak and Sreehari 8 introduced the 
BGD with pmf g(.,.) of a random vector 
(X ly Xi) given by 


20 
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g (m,n) 


f m + n^ 


V n J 


P? Pi Po ,*n.n= 0,1, 2, ..,(13) 


0 < p^Pi <1,0 <p 0 =\-p\ -p 2 < 1. 


where 


A(Pi,P 2 ) 


ZZ a m, n Pi m p n 2>PhP2,a m . n 


m~0n=0 


> 0 . 


The marginal distributions are geo- 
metric with EX i =p\/p 0 and EX 2 = p 2 / po. 

We give an illustration of this BGD in 
Statistical Quality Control. 

Example : 

In the three-class attributes inspection, a 
unit which meets all the specifications 
(important as well as less important) laid 
down by a quality control engineer is called 
a good unit. A unit is called a bad unit if it 
fails to meet one or more important 
specifications. A unit which is neither good 
nor bad is called marginal. Let p 0 , p\ and p 2 
(=1- p 0 - P\) be the probabilities of a unit 
inspected being bad, marginal and good 
unit respectively. Suppose that during the 
unit by unit inspection process, the number 
X\ of marginal units and the number X 2 of 
good units that pass the inspection before 
the first bad unit is encountered are 
recorded. Then the probability mass 
function of (XuX 2 ) is given by (13). 

Phatak and Sreehari proved two 
characterization results concerning the BGD at 
(13) based on properties 1 and 2. 

Definition 7 : A random vector (X\, Xi) 
has a bivariate power series distribution if 
its pmf is given by 

n n m n n 

f(m,n) = , m,«=0,l,2,... (14) 

AP\>Pi) 


Theorem 6 : Let (X\, Xi) have the 
bivariate power series distribution given at 
(14). Then (X u X 2 ) has BGD(P & S) given 
at (13) if and only if 

E(X i )=p i /(l-p l -p 2 ),i=l,2. (15) 

Remark : This result is an extension of 
the property 1 of the geometric distribution. 
The next result is an extension of the 
convolution property 2 of geometric 
distribution. 

Theorem 7. Let (X\ Xi) and {Y\,Yi) be 
independent identically distributed random 
vectors with pmf g(.,.). Then g(.,.) satisfies 
the relation at (13) if and only if 

P [X\+Y\ = m, X 2 +Y 2 =n] = C (m+n+l)g(m,n) 

' (16) 

for all m, n= 0, 1 , 2 ... where C is a constant. 

Phatak and Sreehari proved the 
Theorem 7 by the method of induction. 
They mentioned that (16) is equivalent to a 
partial differential equation satisfied by the 
probability generating function correspond- 
ing to g(.,.) and that by solving that partial 
differential equation one might get a 
simpler proof which can be extended to 
obtain a characterization of multivariate 
geometric distribution given below. 

Definition 8 : A random vector X = 
(X\,...Xi) has a multivariate geometric 
distribution (MGD) if 
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P{X J = n j! j=l,...,k) = 



k 

where 0< p s < 1 , j= 0, 1 ,2 . . . ,k, £ = 1 • 

7=0 

Srivastava and Bagchi 13 proved a 
characterization result concerning the 
above MGD based on the convolution 
property (extension of Theorem 7) solving 
a partial differential equation satisfied by 
the probability generating function of the 
MGD as conjectured by Phatak and 
Sreehari. 

Phatak and Sreehari considered certain 
bivariate versions of the LMP at 1 and the 
property 3 of the univariate geometric 
distribution. 

Extension of LMP at 1 to two 
dimensions 

One may consider 
P[X i = m+l,X 2 =n+l \Xi>l,X 2 >l] = 
P[X x =m,X 2 =n],m,n= 0,1,2... (17) 

as an extension of the LMP at 1 . 

The BGD(P & S) does not satisfy (17). 
Nagaraja 23 proposed a different extension 
of LMP at 1 to the bivariate case which 
indeed characterizes the BGD(P & S). 
Nagaraja proved the following : 

Theorem 8. A random vector (Aj, X 2 ) 
has BGD(P & S) if and only if its pmf 
g satisfies the property 

22 


g ( m,n ) =c x g{m-\,n ) + c 2 g(m, n-1), m,n 
= 0,1,2,..., (0,0) (18) 

and g{m,n ) =0 if m< 0 or n< 0, 

where 0 <c x ,c 2 < c, +c 2 <1. In this case 
g(0,0)=l-c,-c 2 . 

Nagaraja also proved the following to 
strengthen a result of Phatak and Sreehari. 

Theorem 9. For a random vector (X\ JC 2 ) 
any two of the following statements 
together imply the remaining.. 

(a) Condition (17) holds 

(b) Xj and X 2 are independent. 

(c) X\ and X 2 are geometric r.v.s with P(X \ 
=0,X 2 =0) = P(X=0) P(X 2 =0). 

Remark : It may be mentioned that 
Result 3 of Nagaraja in this context is 
incorrect. 

Srivastava and Bagchi 13 obtained a 
different characterization result for the 
BGD(P & S ) using a bivariate power series 
distribution and a random version of 
condition (17). Their result is a 
generalization of a result of Zijlstra 24 , and 
makes use of the previous theorem of 
Nagaraja. 

Theorem 10. Let X= {X\, X 2 ) and Y = 
(Y\, Y 2 ) be independent random vectors and 
let Y have the bivariate power series 
distribution given at (14). Then AT is BGD(P 
& S ) if and only if 

P [X= Y, Y* 0 ] = c, P[Xi= Y\ -1,X 2 =Y 2 ] 
+c 2 P[X 1 =Y u X 2 =Y 2 - 1] 

for all p\, p 2 such that 0< p\, p 2 < M where 
M, c\, c 2 are positive and Cj+ c 2 <1. 
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Phatak and Sreehari considered the 
following extension of Property 3 of 
univariate geometric distribution due to 
Puri 25 discussed earlier. 

For some random vector Y 

P(Y< n)-P(X + Y<ri) =J3P(X+ Y=n) 

\ 

holds for all n = ( n x ,n 2 ), «,,n 2 = 0,1,2,... 
where p is a constant. 

It turns out that there is no Y which 
satisfies the above condition when X has 
BGD(P & S ). Nagaraja considered an 
alternative bivariate extension of the 
condition at (6) and proved the following. 

Theorem 1 1 : A random vector X = (X\, 
Xi) has BGD(P & S ) if and only if for 
some random vector Y = (Y\, Yi) 


P(X\=m) = Po 


1 -Pi 


Pi 

J-P 2 


m- 0,1,2... 


where 0<p o , p\,p 2 < \,p 0 + P\ + p» 2 = l- Then 
the joint distribution of (X\, Xi) is given by 
(13). That is, (Xi , Xi) is BGD(P & S). 

The proof appears to be incomplete and 
as such the problem considered is open. 

A similar result in the case of BGD(N 
& N) will be discussed later. 

Bivariate geometric distribution BGD 
(N & N). 

Nair and Nair 9 considered the following 
extension of mean residual life function p 
(t) given in definition 3 . 


P (m,n) = (pi (m \ n), p 2 (« I m)) 


P(Y= n) = (1+p) P (X + F = /i ) - ep P(X + 
Y = B -70- (1-0) p P(X+Y = n-Ii) 


where I\= (1,0), h= (0,1), 0 and p satisfy 
Po = ( 1 + P )'' 5 A =0 P Po 311(1 P 2 = d- 0 )P Po- 

Nair and Nair 20 proved interesting 
results similar to Theorem 3 for the BGD(P 
& S) and BGD(N & N). 

Theorem 12. Let X\ and Xj be two r.v.s 
for which the conditional distribution of X\ 
given X 2 =n is given by 


P(X x = m\x 2 =n) = 


fm + n' 


v m J 


Pi (1 -Pi) n+ \ 
m= 0 , 1 , 2 ,... 


and the marginal distribution of X\ is given 
by 


\Li(m\n) = E[Xi -m\Xi>m,X 2 >n] (19) 

and 

H 2 (n\m) = E[X 2 -n\ Xi>m,X 2 >n] ( 20 ) 
m,n = 0,1,2,... 

Recall that the property that p(n) is a 
constant for all n- 0 , 1 , 2 ..., characterizes a 
geometric distribution. Nair and Nair 
proposed that if pi (m | n) and p 2 (n | m) are 
independent of m and n then the joint 
distribution of (X\, Xi) may turn out to be a 
bivariate version of geometric distribution. 
But it turns out that X\ and X 2 are 
independent geometric r.v.s. They modified 
the constancy of pi (m \ n) conditions which 
led to the BGD(N&N). Suppose 
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p\ (m | n) = a x (n) is independent of m (21) 

and 

p 2 (« | m) = a 2 (m ) is independent of n. (22) 

They proved the following result. 

Theorem 1 3 : A random vector (X \ , X 2 ) 
has the BGD given by the joint survival 
function 

R(m,n) = p? p\ 6 mn , m,n= 0,1,2.. . (23) 

where ti<p\<pz<\, O<0 <1, 1-0 < (1 -pfi ) (1- 
p 2 Q ) if and only if pi (m | n) and \x 2 (n \ m) 
satisfy (21) and (22) where a\(n) and a 2 (m) 
are non-increasing and ai(0) =p, / (1 -pi), 

/= 1 , 2 . 

In view of the equality of conditions (4) 
and (5) for a geometric r.v. one looks for a 
bivariate version of the LMP which might 
also characterize the BGD given at (23). 
Such a result is proved by Nair and Nair 9 in 
the following theorem. 

Theorem 14 : A random vector (X\, Xz) 
has the BGD given at (23) if and only if the 
following conditional lack of memory 
properties 

P[ X x > m+k | Xi >m, X 2 >n] = P[X x >k \ X 2 >n\(2A) 
and 

P[X 2 > n+k | Xi >m, X 2 >n] = P[X 2 >k \ X x >m]{25) 

€ . 

hold for all non-negative integers m,n and k. 

Analogous to Theorem 3 and Theorem 
12 we have the following distributional 
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property characterizing BGD(N & N) 
proved by Nair and Nair 20 . 

Theorem 15 : Suppose 
P(Xi=m\X 2 =n)= 

P :'e-' a - Pl r [o-p, e- +, )a-p a 0’ +l ) 

+ 0-1], m = 0,1,2,,. . 

and 

P(X\=m) = (l-pi) p x , m=0,l,2,... 
where 

0<p,,p 2 <1,0 <©<1,1-9 <(l-p, 0)0 -A0). 

Then (Xu X 2 ) has the BGD given at 
(23). 

Nair and Asha 10 discussed the IFR (and 
DFR) properties of BGD(H), BGD(N & N) 
and BGD(N& A). 

Bivariate geometric distribution, BGD 
(A & S) 

We now discuss a reliability problem 
which leads to the BGD considered by 
Asha and Sankaran. 

Example : 

Let X\ and X 2 be the r.v.s representing 
the lifetimes (measured in the number of 
cycles of survival) of two components A 
and B in a two component system. Suppose 
that the two components are subject to 
shocks Si and S 2 respectively. The shocks 
occur independent of each other at most 
once in a cycle. Suppose Z\ and Z 2 are the 
numbers of the cycles in which shocks Si 
and S 2 occur causing failure of components 
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A and B respectively. Clearly component A 
fails first in cycle m if Z\=m and Z 2 >Z\ and 
similar is the case with component B. After 
the failure of either component, the system 
will continue to operate till a further shock 
(to the other component) occurs. Let 
Z\ (and Z*) be the number of the cycles in 
which shock Si (S 2 ) occurs causing failure 
of component A (and B) if Z\<Z 2 (Z\>Z 2 ). 
Then the component lifetimes are 


bivariate LMP considered at (9). Further, 
the BLMP, under certain additional 
assumptions concerning the marginal 
distributions, characterizes the above BGD. 

The marginal distribution of Xj is a 
mixture of two geometric distributions and 
is given by 

P(X r m) = 


X = 



if Z, < Z 2 
ifZ x >Z 2 


and 


f z 2 </z,>z 2 

[Z, + Z 2 if Z, < Z 2 . 


P f\- p'Mr 


-(piPiY 


Pj ~ P\Pi 

Pjf_ - A-yX 1 -P))~ iP) ~ Pi AX 1 ~ Pj) 
Pj-PxPi 


(27) 


m = 0,1,2,... 


Under the assumption that Zi, Z 2 , 
Z\ and Z* are independent geometric r.v.s 

with expectations p\, p 2 , pi and p\ , the 
joint distribution of (X\, X 2 ) is given by 

P(X x =m,X 2 =n) = 


They obtained a characterization result 
for the BGD given at (26) through bivariate 
mean residual functions defined as follows: 

Definition 9 : For two r.v.s X\ and X 2 
the bivariate mean residual life is the vector 
(9 (t), dfn | m), Ofm \ n )) where 


(A Pi )" (Pi Y ” -1 Pi (1 - Pi )(1 - Pi ) 

if m>n 

{PiPiY (p'iY~ m ~ X Pi(y~ P\)(}~ Pi) 

if m<n 

(PiPiY 0-A)0“A) if m=n 

(26) 


0 (t) = E[ minCJSTi, X 2 ) - t\x x >t, X 2 >t] 

Qi(n | m) = E[X\-n | X\>m, X 2 =n\ for m > n 
02 (m | n ) = E[X 2 -m | X\-m, X 2 >n] for n>m. 


m,n = 0 , 1 , 2 ,.. 
where 

Pi >Pi * P 1 P 2 ,0<p l ,p 2 ,pl,p' 2 <\. 

Asha and Sankaran observed that 
min^i^) is geometrically distributed with 
mean pip 2 . The distribution satisfies the 


Since failure rate and mean residual life 
times are related it is natural to expect that 
9 (t), 6\ (n | m) and 9 2 (m | n) are related to 
r(t), r\(m \ n) and r 2 (n | m) defined by Sun 
and Basu 17 . Asha and Sankaran established 
the relationship. Further, they proved the 
following characterization result concern- 
ing the BGD given at (26). 
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Theorem 16. Let (XxJC-T) be a random 
vector with marginals given at (27). If 
P[Xx>X 2 \mm{Xi > Xi)=z] and P[X i< X 2 \ 
are constants, then the 
bivariate mean residual life (0 (t), 6\(n | m), 
& 2 .(m | n)) is a constant vector if and only if 
(Xi^f 2 ) has the BGD given at (26). 

Concluding Remarks 

There are several more characterizing 
properties of univariate geometric distribu- 
tion. It will be interesting to formulate their 
analogues in two dimensions and examine 
whether they lead to any of the bivariate 
geometric distributions considered here or 
lead to some new distributions which may 
have a useful role in stochastic modelling. 
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Abstract 

Fundamentals and recent advances in electro- 
chemical Supercapacitors, which are recently 
emerged as an alternative power sources, have 
been discussed. A general discussion on different 
types of Supercapacitors namely: electrical double 
layer capacitors (EDLCs) and electrochemical 
redox Supcrcapacitors has been presented in view 
of their configuration, materials and charge 
storage mechanisms. Particular attention and 
emphasis are given to the conducting polymers 
based redox capacitors and their different types 
(I, II and III) of the configuration according to 
three different possible combinations of p- and n- 
doped conducting polymers. Recent investigations 
in our laboratory on the development of different 
Supercapacitors have been briefly given. Various 
possible applications of Supercapacitors as a 
power source have also been briefly described. 

(Keywords : Supercapacitor/ Polymer Electrolytes/ 
Gel electrolytes/conducting polymers. 

Introduction 

Capacitors are important electrical 
storage devices, which are useful electronic 
components in high frequency circuit 
applications such as electronic commu- 
nication, telephone exchanger etc. and 
medium frequency circuits like power 
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filters, etc. In general, the existing 
capacitors can be grouped into three 
categories: (a) conventional capacitors 
(ceramic, paper, thin film), (b) electrolytic 
capacitors, and (c) electrochemical capa- 
citors. The configuration of the conven- 
tional capacitors are simple in which a mica 
/paper/ ceramic film is sandwiched between 
two metallic electrodes. The configurations 
of other two types of capacitors are 
described as follows. The characteristics 
and possible applications of different types 
of conventional and electrolytic capacitors 
are summarised in Table 1. 

There exist two types of electrolytic 
capacitors, which are commonly available 
in market namely: 

(i) Aluminium electrolytic capacitors, and 

(ii) Tantalum electrolytic capacitors 

These capacitors use liquid as well as solid 
electrolytes. Due to various advantages of 
solid electrolytes over the liquid electrolytes 
such as no leakage, improvement in 
temperature and frequency response, heat 
resistance and possibility of sleek design for 
miniaturisation, the solid electrolytes based 
electrolytic capacitors are more popular. The 
sketch diagrams of tantalum and aluminium 
solid electrolyte capacitors are shown in Fig.l 
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Table 1- Specification and possible applications of conventional and electrolytic capacitors 


Capacitor types 

Specification 

Application 

Paper 

-400 V, <0.1 pF 

High frequency electronic circuits 

Mica 

-100 V, <1 pF 

High frequency electronic circuits 

Ceramic 

-50 V, <0.1 pF 

High frequency electronic circuits 

Thin film 

-50 V, <0.01 pF 

High frequency electronic circuits 

Aluminum Electrolytic (Liquid) 

- 50 V, >1 pF 

Medium frequency electronic 
circuits. Power supply filters, etc. 

Tantalum Electrolytic 
(Liquid) 

-50 V,>1 pF 

Medium frequency electronic 
circuits, Power supply filters, etc. 

Aluminum Electrolytic 
(Solid) 

-25 V,>1 pF 

Medium frequency electronic 
circuits, Power supply filters, etc. 

Tantalum Electrolytic 
(Liquid) 

-25 V, >10 pF 

Medium frequency electronic 
circuits, Power supply filters, etc. 


Ta- 


( + ) 


Liquid 

Electrolyte 

or 

Solid MnO-, 


I 


^Counter 

Electrode 

<-) 


Ta 2°5 


(A) Tantalum Electrolytic Capacitor 



(B) Aluminium Electrolyte Capacitor 


Electrochemical Supercapacitors : 

Over the past one decade, a great 
interest in ‘Supercapacitors’, commonly 
referred as ‘Electrochemical Capacitors’ 
has been stimulated worldwide due to the 
growing demand for power sources of 
transient high power density. 1 These are 
potential devices to work as an alternative 
power source for low energy density 
applications e.g. electronic and medical 
equipment, whereas due to their high rate 
capability, they are also considered to be 
excellent candidates for use in high power 
applications particularly for load levelling 
in electrical vehicles. Electrochemical 
Capacitors, are unique energy storage 
devices which have the following specific 
characteristics: 

• Large value of capacitance (-20-2000 
times per gram or per cm 2 of 
electroactive materials). 


Fig. 1— The configuration of tantalum and aluminium 
electrolytic capacitors. 


High power storage/delivering ability. 
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Historically, Helmoltz in 1879 first 
considered the storage of charge at 
electrode-electrolyte interface and regarded 
the double layer as equivalent to a parallel 
plate capacitor. The formation of Helmoltz 
double layer, which acts as a parallel plate 
capacitor is shown schematically in Fig. 2. 
Typical values are in the range of 20-40 pF 
cm' 2 for electrode-liquid electrolyte 
interfaces. Rayleigh (1966) reported a 
substantially higher capacitance of 200 pF 
cm' 2 for the interface of platinum and 
AgBr. This higher value was explained in 
terms of the adsorption of a layer of ions on 
the electrode surface. In 1980, NEC 
Corporation (Japan) reported and 
commercialised the Electrical Double 
Layer Supercapacitors, giving hundreds of 
mF cm* 2 specific capacitance, composed of 
large surface area activated carbon 
electrodes and sulfuric acid electrolyte with 
a brand name ‘SUPERCAPs.’ Thereafter 


the Supercapacitors of various types 
employing different electrodes and 
electrolytes have been reported in literature. 
For details, readers are referred to the 
reviews / books. 1 ' 5 

The configuration of Supercapacitors is 
almost similar to that of batteries, the 
difference being that the former consists of 
electrolytes sandwiched between two 
blocking / active electrodes with translation 
of charge only, whereas a battery employs 
asymmetrical different electrode materials 
(anode and cathode) involved in interfacial 
reactions leading to an overall electro- 
chemical process for the cell. 

In recent developments, two types of 
electrochemical capacitors have been 
investigated which are based on two 
fundamental charge storage mechanisms. 
These are classified as follows: 


Supercapacitors — 


Electrical Double 
Layer Capacitors 
(EDLCs) 


Redox 

Supercapacitors 


I — Electrostatic Capacitance 


I — Electrodes : Carbon 


Earadic Pseudo-capacitance 

Insertion type 


Electrodes 


compound (e.g.Ru02 
NiO etc.) 


Conducting Polymers 
(e.g. polypyrrole, etc.) 
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®/+ Cation - Anion 


Fig. 2- General representation of Helmoltz double 
layer and its potential profile. 

Electrical Double Layer Capacitors 
(EDLCs) : 

The capacitors, which employ the 
carbon or other similar materials as block- 
ing electrodes are referred as “electrical 
double layer capacitors (EDLCs).” 
Electrodes, which do not participate in any 
electrochemical reaction including charge 
transfer at . the electrode-electrolyte 
interface (within the voltage limit of 
electrochemical stability range of elec- 
trolyte), are able to accumulate (store) 
charges under the application of electric 
potential. These electrodes are regarded as 
polarisable or blocking. In general, large 
surface area carbon electrodes are being 
used in EDLCs. Several electric double 
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layer capacitors using large surface area 
electrodes with different types of carbons 
have been developed with liquid and 
polymer / gel electrolytes. 6 ' 11 In particular, 
activated carbon powders / fibres / fabrics, 
which can be produced with large surface 
area up to 2000-2500 m 2 g' 1 , have attracted 
much attention for the fabrication of double 
layer capacitors due to their corresponding 
very high capacitance values of 400-500 
Eg' 1 . The charge storage mechanism in 
EDLCs is electrostatic resulting from the 
separation of charged species at the 
interfacial double layers. The charge 
storage mechanism and the operational 
principle of EDLCs are shown schemati- 
cally in Fig. 3. 

The capacitance, C and the charge, Q, 
accumulated in the electrical double layer 
formed at the interface between the 
blocking (polarisable) electrodes and 
electrolyte is expressed as 4 : 

Capacitance C = [jsMrorjiiS' (1) 

Accumulated charge 

Q- [js/(47ca)nS'](2\|/ 1 ) (2) 

where s is dielectric constant of the 
electrolyte, a is distance from the electrode 
interface to the centre of ions and S is 
electrode surface area. 

Redox Supercapacitors 

Capacitors, which employ the electro- 
active materials such as insertion type 
compounds (noble metal oxides e.g. R.UO 2 , 
NiO etc.) 12 ' 14 or conducting polymers (e.g. 
polypyrrole, polythiophene etc.) 15 ' 25 as 
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Fig. 3- Schematic structure of electrical double layer capacitor and its potential profile. 


active electrodes are referred as ‘redox 
supercapacitors. ’ In redox supercapacitors a 
fast faradic charge transfer process takes 
place which gives rise to a ‘pseudoca- 
pacitance. ’ Recently a new class of supra- 
molecular materials (e.g. 1,5 Diamino- 
anthraquinone oligomer, etc.) have been 
introduced as promising electrodes for 
redox supercapacitors, which exhibit high 
specific capacity with aquous and non- 
aquous liquid electrolytes. 26, 27 

Particularly, conducting polymers have 
been shown to be potential electrode 
materials for redox supercapacitors because 
they offer several advantageous properties 
over the other noble metal oxide/carbon 
electrode materials. These include their low 
cost, easy chemical or electrochemical 
preparation in desired extended area thin 
films and large capability of charge storage 
SCIENCE LETTERS, VOL. 27, NO. 1 & 2, 2004 


(-500 Cg' 1 ) throughout their volume. 
Particularly for electrochemical capacitors, 
the ability to charge (or discharge) through 
the whole volume is most attractive aspect 
of conducting polymers compared to the 
large area carbon materials. It is well 
known that conducting polymers exist in p- 
doped (oxidised) and n-doped (reduced) 
states. In view of this fact, Rudge et al 
proposed three schemes of configuration of 
conducting polymers based redox 
capacitors, classified namely as types I, II 
and III. Two identical p-dopable polymers 
are employed as active electrodes in type I 
capacitors, whereas type II capacitors are 
comprised of two different p-dopable 
conducting polymers sandwiching the 
electrolyte materials. In type III capacitors, 
a configuration of n- and p- doped forms of 
the same polymers can be used. The 
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classification of different types of redox 
capacitors under these three categories is 
based on the difference in potential range in 
which the conducting polymers become 
doped or dedoped, resulting in an increase 
of voltage of the charged capacitors. The 
type 1 capacitors are typically limited to ~ 1 
volt whereas the type III can be charged to 
~3V. 

Electrolytes 

In most of the studies on the double 
layer or redox supercapacitors, liquid 
electrolytes have been used as separators. 
These include acetonitrile (ACN)- 
(CH 3 ) 4 NBF 4 , ACN-(C 2 H 5 ) 4 NBF 4 , ACN- 
(C 4 H 7 ) 4 NPF 6 , ACN-LiC10 4 , propylene 
carbonate PC-LiC10 4 , PC-(C 2 H 5 ) 4 NBF 4 etc. 
Liquid electrolytes based supercapacitors 
are similar to liquid electrolytes based 
batteries in that they are associated with the 
well known disadvantages of corrosion, self 
discharge, low energy density and bulky 
design. A limited number of solid-state 
capacitors have been studied with polymer / 
gel electrolytes e.g. poly (ethylene oxide) 
(PEO)-PC-(C 2 H 5 ) 4 NBF 4 , Poly (ethylene 
glycol) (PEG)- (C 2 H 5 ) 4 NBF 4 , poly (methyl 
methacrylate) (PMMA)-PC-ethylene car- 
bonate (EC)-LiC10 4 , polyvinyl alcohol 
(PVA)-H 3 P0 4 , PEO - PEG - LiCF 3 S0 3 etc. 
All solid-state supercapacitors are in the 
early stage of. development and require 
improvements in their capacitance values, 
reduction of internal resistance and 
enhancement in the electrochemical sta- 
bility range of potential. Solid electrolytes 
(polymer / gel electrolytes) as separators 
possess some limitations in terms of lower 
conductivity, which leads to the higher 
internal resistance and hence lower power 
density. The electrochemical potential 
ranges of solid electrolytes are also 
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relatively lower as compared to particularly 
organic liquid electrolytes, which results 
the lower energy density of the capacitors. 
Some important examples of electrical 
double layer and redox supercapacitors 
which employ different liquid / polymer / 
gel electrolytes are listed in Table 2. 

Physical Techniques for Supercapacitor’s 
Characterisation 

Most of the experimental techniques, 
which are commonly employed for the 
charcterisation of electrochemical cells 
(including batteries, sensors, electro- 
chromic devices, etc.) are equally suitable 
for the supercapacitors. These include a. c. 
impedance spectroscopy, linear sweep 
voltammetry, charge-discharge methods, 
self discharge, prolonged cyclic tests etc. 
The performance parameters associated 
with the supercapacitors can generally be 
obtained by discharging a fully charged 
device through a load resistance. The 
parameters along with their expressions and 
units are listed as in Table 3. 

A. C. Impedance Spectroscopy 

An electrochemical capacitor cell, as 
described previously, is generally 
fabricated from two identical electrode 
systems with a liquid or a solid electrolyte 
separating them. Impedance spectroscopy is 
an attractive method for characterising the 
electrochemical behaviour of such cells 
because of the direct connection between 
the real system and its idealised equivalent 
circuit consisting of discrete electrical 
components (R, C and L) in their series and 
parallel combinations. As described earlier, 
the real systems in the present context are 
electrochemical capacitor cell which 
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employ either blocking / polarisable planar geometry (e.g. high density graphite 

electrodes (called double layer capacitors) sheets, glassy carbon sheets, etc.) or large 

or electroactive electrode materials (called surface area porous electrodes (activated 

redox supercapacitors). The blocking carbon powder/ fabrics/ fibre etc.), 

electrodes used in the capacitors have either 

Table 2- Characteristics of some important EDLCs and redox supercapacitors 

Capacitor Ceils 

Capacitance (C) 
(Fg 1 ) (mF cm' 2 ) 

Energy 

Density 

(Wh Kg') 

Power 
Density 
(kW Kg') 

Ref. 

Electrical Double Laver Capacitors 

ACF | Nylon 6-10 + 2 H 3 P0 4 . | ACF 

70 


1.4 


7 

Ni-ACF | PC-LiC10 4 | Ni-ACF 

95 


- 

- 

8 

ACF | PMMA-PEO-(C 2 H 5 ) 4 NBF 4 - PC 1 ACF 

20 

570 

16.2 


6a 

ACF 1 PAN - (C 2 H 5 ) 4 NBF 4 - PC 1 ACF 

16 

370 

- 

- 

6b 

ACF 1 PE0+PEG+LiCF 3 S0 3 1 ACF 

4 

20 

- 

- 

10 

ACF 1 PVA + H 3 P0 4 | ACF 

70-90 

360-470 

14-18 

- 

10 

AC | PVdF-HFP-EC-PC-(C 2 H 5 ) 4 NBF 4 |AC 

. 123 

- 

- 


9 

HD-ACF | Aqueous KC1 1 HD-ACF 

28 

- 


- 

ii 

Redox Supercapacitors 

Ru0 2 .xH 2 0 |H 2 S0 4 (aq.) | Ru0 2 .xH 2 0 

720 

- 

26.7 


12 

Ni-NiO | 1 M KOH (aq.)| Ni-NiO 

50-64 

- 

7.2-11.0 

■ 4.0-17.0 

13 

PPy 1 PMMA-EC-PC-LiClO„ 1 PANI 

13 

- 

4.0 

0.7 

21 

PPy | PVA + H 3 P0 4 | PPy 

39-84 

1. 5-5.0 

12.0 

- 

23 

PMeT | PE0-PEG-LiCF 3 S0 3 1 PPy 

25-30 

3. 0-3. 5 

10.0-12.0 

- 

22 

pDTT | PEO-PC-(C 2 H 5 ) 4 NBF 4 | pDTT 

25-35 

5.0-7.0 

3. 0-7.0 

0.2 

19 

PEDOT | ACN-LiC10 4 1 PEDOT 

- ■ 

- 

1. 0-4.0 

2.5 

20 

PPy / graphite | KC1 (aq.) | PPy/ graphite 

400 

- 

- 

- 

28 

Li-PANI | ACN-(C 2 H 5 ) 4 NBF 4 1 Li-PANI 

60-100 

- 


- 

29 

pDAAQ | PC-(C 2 H 3 ) 4 NC10 4 1 pDAAQ 

- 

- 

25.0-46.0 

10.2-30.5 

26 

pDAAQl PMMA-EC-PC-(C 2 H 5 ) 4 NC10 4 | pDAAQ 

3. 7-5. 4 

125-184 

92-135 


30 

ACF : Activated Carbon Fabrics, Ni-ACF : Nickel 

loaded ACF, HD-ACF 

: High Density ACF, AC 

: Activated 


Carbon Powder, PPy : Polypyrrole, PAN1 : Polyaniline, Li-PANI : Lithium doped PANI, PMeT : Poly (3-methyl 
thiophene), pDTT : Poly (thieno [3,4-b : 3’,4’-d] thiophene, PEDOT : Poly (3,4 ethylene -dioxythiophene), pDAAQ 
: poly (1,5-diaminoanthraquinone). 
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Table 3- Performance parameters of Supercapacitors. 


Parameters Expression 

Unit 


Capacitance 
Specific Capacitance 


C= q/V 

C’ = C/M 
C” = C/A 


Farad ‘F’ 

Fg' 


F cm 


-2 


Stored Charge 


Energy Density 


Power Density 



u J 


2 \ 


P = 


M 
_i 
2M 

r 2 

mr l 


v_ 

R,.. 


dt 


-CV 2 


A schematic representation of the 
impedance behaviour of an ideal capacitor 
cell is illustrated in Fig. 4(a). The 
impedance pattern shows a steep rising 
dispersion line overlapping with imaginary 
axis of the impedance plot. The impedance 
pattern of a real capacitor (Fig.4 b) shows a 
perfect semicircle representing a parallel 
combination of a bulk resistance Rb 
(resistive component owing to the ion 
migration owing to the bulk of electrolytes) 
and geometrical capacitance C g (due to 
dielectric polarisation of electrolyte under 
the influence of an electric field). The 
interfacial polarisation leads to the 
formation of double layer capacitor at the 
two identical interfaces, which is 
represented by a steeply rising spike in the 
low frequency region of the impedance 
spectra (Fig. 4 b). 

The equivalent circuit of capacitor cells 
with either rough or porous electrode in 
contact with the electrolytes can be 
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Coulomb 


Jg 1 or Wh Kg' 1 


Jg's 1 or W Kg' 1 


represented in the form of a transmission 
line type distributed circuit which gives 
variable time constant. Conducting polymer 
electrodes (e.g. polypyrrole, polytiophene 
and its derivatives etc.) also show a similar 
pattern of impedance behaviour as the 
porous electrodes. The region of the 
impedance spectrum associated with high 
capacitance at low frequencies is often 
interpreted in terms of a distributed R and 
C network due to ionic penetration 
throughout the polymer films or the porous 
electrodes. 

The overall capacitance of the cell can 
be evaluated from the relationship : 


where 2” is the imaginary part of the total 
impedance measured in low frequency range 
(100 mHz- 1 mHz). 
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Fig. 4- Schematic representation of complex impedance Fig. 5- Schematic representation of linear sweep 

plots for (a) ideal capacitor and (b) real capacitor. voltammograms for (a) ideal capacitor and (b) 

real capacitor. 


Linear Sweep Voltammetry 

The linear sweep voltammetry is 
another important experimental probe to 
characterise the capacitors. The cyclic 
voltammogram of an ideal capacitor 
(having negligible internal resistance) 
shows a perfect rectangular shape (Fig. 5a) 
and it is scan rate dependent. The deviation 
from the rectangular shape is observed 
(Fig. 5b) for real capacitors which is 
represented as a series combination of 
internal resistance R and ideal capacitance 
C. The value of capacitance can be 
evaluated from the expression : 

(4) 

s 
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where ‘f is the voltammetric current and ‘4’’ 
is scan rate. 

Charge-Discharge Methods 

Information about an electrochemical 
system is generally gained by applying an 
electrical perturbation to the system and by 
observing the resulting changes in the 
characteristic behaviour of the system. In 
impedance spectroscopic measurements an 
a.c. perturbation is applied, whereas a 
voltage ramp is applied in the case of linear 
sweep voltammetry. These methods 
have already been described earlier. In 
the charge-discharge technique for 
characterisation of a capacitor, it is 
necessary to consider the behaviour of a 
real capacitor system represented by an 
ideal capacitor C e connected with the 

35 





LEAD ARTICLES/OVERVIEWS OF NEW DEVELOPMENTS 


internal resistance Rj in series. The 
perturbation can be either by potential step 
(called the potentiostatic method) or by 
current step (called the galvanostatic 
method). 

(a) Potentiostatic method (Voltage step) 

When a potential step is applied across 
the capacitor cell, the capacitor charges and 
an exponentially decaying current is 
obtained. The behaviour of the current ‘i’ 
with respect to time, on application of the 
potential step Vo, can be expressed as : 



where ‘x (= RC )’ is the time constant. 
The capacitor cells start discharging when 
the operation is reversed and the pattern of 
discharge curve is similar to the charging 
curve but. with the reverse polarity of the 
current. 

The charge-discharge curves can be 
integrated numerically to obtain the total 
charge stored or released from the cell. The 
capacitance ‘C’ can be evaluated by 
dividing the total charge by the applied 
voltage Vo i.e. 


(b) Galvanostatic method (Current step) 

When a capacitor cell is charged by a 
constant current the total voltage ‘P 
developed across the cell can be expressed 
as : 


As the current is constant : 


V = i 


R, 


+ 



( 8 ) 


Thus, for a current step, the voltage across 
an ideal capacitor cell, whose capacitance 
value, C e , is expected to be constant, 
increases linearly with time. A schematic V- 
t curve showing the linear variation is shown 
in Fig. 6(a). After the full charging of the 
capacitor at the maximum permissible or 
desired voltage level, it can be discharged 
galvanostatically. The voltage decrease 
across the cell is also linear with respect to 
time (Fig. 6a). In the real capacitors, The 
initial sharp jump in voltage while charging 
or discharging of the capacitor is observed 
(Fig. 6b), which is owing to an ohmic (iR) 
loss due to internal resistance of the cell. 



V = iR, 
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(7) 


Fig. 6- Schematic representation of galvanostatic 
charge-discharge curves for (a) ideal capacitor 
and (b) real capacitor. 
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The capacitance of the cell can be 
evaluated from the linear region of the 
curve following the relation : 


where AV is the change in voltage for At 
interval of time. The charge-discharge 
experiment can be carried out for several 
cycles (10 3 to 10 5 cycles) to test the cyclic 
efficiency of the cell. It may be noted that 
the cyclic efficiency of supercapacitors is 
100-1000 times larger than that of the 
rechargeable (e.g. Ni-Cd, Ni-metal hydride 
etc.) batteries. This aspect makes the 
supercapacitors more attractive as power 
source over the rechargeable batteries. 

The 'coulombic efficiency can also be 
calculated easily from the results of 
galvanostatic charge-discharge experiment. 
When constant current is used in charging 
and discharging the coulombic efficiency 
‘rj’ can be described by : 


n= — — jci oo% (10) 

At r 


where Ate and A to are time intervals for 
charging and discharging, respectively. 

Some Experimental Studies 

Some important results recently 
obtained in our laboratory on different 
types of double layer and redox 
supercapacitors have been discussed. 
Particular emphasis on all solid state 
supercapacitors based on polymer / gel 
electrolytes have been given. 
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Polymer / Gel Electrolytes 

Various polymer / gel electrolytes have 
been synthesised and used to fabricate 
different types of solid-state supercapaci- 
tors. These include polymer electrolytes, 
PVA-H 3 PO 4 aqueous blend, PEOi 5 .NaClC >4 
+ PEG and polymeric gel electrolytes, 
PMMA-EC-PC-NaC10 4 and Polyvinyli- 
dene fluoride-co-hexafluoropropylene 
(PVdF-HFP)-EC-PC-alkyl ammonium 
salts. Films (-200-300 pm thick) of the 
polymer/gel electrolytes were prepared 
using standard ‘Solution Cast’ technique. 
Details of the preparation methods of these 
films have been described by us 
elsewhere. 24, 25 

Fig. 7 shows the typical electrical 
conductivity profile of polymer and gel 
electrolytes with respect to temperature for 
the temperature range 0°C to 150°C. The 
polymer electrolytes based on PVA and 
PEO are generally obtained in the form of 
free standing thin films. The polymeric gel 
electrolytes based on PMMA and PVdF- 
HFP, used in the present work, have also 
been observed in the form of free standing 
mechanically stable thin films. The 
electrolytes show the conductivity of the 
order of 10’ 4 -10' 3 S cm' 1 at room 
temperature (25°C). This order of 
conductivity of the polymer/gel electrolytes 
is acceptable because it offers low 
resistance (say 100 ohm cm 2 ) when the 
electrolytes are used in the form of thin 
film of smaller thickness (say -200 pm). 

The electrochemical potential window 
i.e. the working voltage range, which is the 
measure of electrochemical stability of the 
electrolytes, is another important parameter 
to be evaluated, which limits the working 
voltage of the capacitors. The 
electrochemical potential window of the 
aqueous based polymer electrolytes like 
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PVA-H3PO4 blend is limited to 1.2 V, 
whereas the PEO, PMMA and PVdF based 
organic polymer / gel electrolytes, used in 
the present studies, offer the potential range 
more than 2.5 V. 



Fig. 7 - Temperature dependence of the electrical 
conductivity of (a) PEO J5 . NaCl0 4 + 50 wt% 
PEG polymer electrolyte and (b) PMMA- 
EC/PC-NaC10 4 gel electrolyte. 

All-Solid-State Electrical Double Layer 
Supercapacitors 

The electrical double layer capacitors 
have been constructed from large area 
activated carbon fabric (ACF, large surface 
area - 2500 m 2 g' 1 , Spectracorp, USA) 
electrodes with the polymeric gel 
electrolytes. The gel electrolytes were 
based on (PVdF-HFP)-EC-PC-tetraethyl 
ammonium tetra fluoroborate (TEABF4). 
Following three different model capacitor 
cells have been constructed using ACF 
electrodes and gel electrolytes of different 
compositions : 

Cell A : ACF | PVdF-HFP (30 wt %) - 
EC:PC(1:1 V/V)-1.0MTEABF 4 | ACF 
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Cell B : ACF | PVdF-HFP (25 wt %) - 
EC:PC (1:1 V/V)-1.0 M TEABF4 | ACF 

Cell C : ACF | PVdF-HFP (20 wt %) - 
EC:PC(1:1 V/V)- 1 .0 M TEABF4 | ACF 

The performance characteristics of the 
capacitors have been characterised using 
a.c. impedance spectroscopy, linear sweep 
voltammetry and prolonged cyclic tests. 

Impedance Studies 

Fig. 8 shows a typical complex- 
impedance spectrum of a capacitor cell C. 
As discussed earlier that the behaviour of 
an ideal capacitor is shown as a straight line 
parallel to the imaginary axis of the 
impedance plots. The impedance response 
of all the cells, in the present studies, shows 
a small semicircular spur in the high 
frequency region followed by steep rising 
portion in the lower frequency region up to 
10 mHz. This steep rising behaviour 
indicates the capacitive character of the 
cells at lower frequency range, whereas the 
high frequency semicircle shows the bulk 
properties of the electrolytes and interfacial 
charge transfer processes. Various 
electrical parameters i.e.. charge transfer 
resistance R ct , the total resistance R and 
capacitance values (evaluated using eqn. 3) 
at frequencies 100 mHz and 10 mHz for all 
the cells A-C are listed in Table 4. A 
comparison indicates that ACF shows a 
highest capacitance value of -28.0 mF cm’ 2 
(equivalent to single electrode specific 
capacitance - 4 Fg' 1 ) at 10 mHz with the 
PVdF-HFP (20 wt. %) - EC:PC (1:1 V/V)- 
1.0 M TEABF4 gel electrolytes. The gel 
electrolyte containing 20 wt.% of PVdF- 
HFP also offers the least value of charge- 
transfer resistance R ct of about 5.0 ohm 
cm 2 . The high value of C and low value of 
R c t is the combined effect of large surface 
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Table 4- Electrical parameters of different EDLCs evaluated from impedance analysis. 


Capacitor 

Ceils 



1 00m Hz 



10mHz 


Ret 

(ohm cm 2 ) 

R 

(ohm cm 2 ) 

C 

(mF cm' 2 )* 

(F g'T 

R 

(ohm cm 2 ) 

C 

(mF cm' 2 )* 

(FgT 

A 

88 

455 

2.3 

0.34 

1800 

4.4 

0.7 

B 

15 

170 

5.0 

0.74 

1170 

7.5 

1.1 

C 

5 

61 

16.4 

2.4 

325 

28.0 

4.1 


* Overall capacitance of the cells 
** Single Electrode Capacitance 

area (~ 2500 m 2 g' 1 ) of ACF electrode and 
highly conducting and mechanically 
flexible gel electrolytes. The capacitance 
value based on 2000 m 2 g" 1 of ACF and 
sulphuric acid liquid electrolyte has been 
reported in literature to be a larger value of 
100 Fg' 1 . 31 The lower value observed in the 
present studies of all-solid state capacitor is 
possibly due to inaccessibility of all of the 
surface area of ACF to form capacitive 
interface with the solid-like gel electrolytes. 

Linear Sweep Voltammetry 

A typical linear sweep voltammograms 
of the capacitor cell C at different scan 
rates between -800 mV and 800 mV is 
shown in Fig. 9. At lower scan rates (say 
10-20 mV sec' 1 ), the votammograms of all 
the cells A-C show an almost featureless 
square shape, which is the basic 
characteristic of a capacitor. The deviation 
from the shape at higher scan rates is owing 
to the equivalent series resistance (ESR), 
which is practically present in series in the 
real capacitors. Further, the response of 
each capacitors cells have been found to be 
scan rate dependent, which is also expected 
for a capacitor cell. 24 



Fig. 8- The typical complex-impedance plot of the 
capacitors cell C: ACF | PVdF-HFP (20 \vt 



Fig. 9- Cyclic voltammogram of the capacitor cell C 
at different scan rates. The scan rates (in mV 
sec' 1 ) are marked on the figures. 
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The values of capacitance have been 
calculated using the eqn. 4 and found to be 
in almost good agreement with the values 
obtained from impedance analysis, 
described above. The cyclic behaviour i.e. 
the stability of the supercapacitors for large 
number of charge-discharge cycles has also 
been evaluated. Fig. 10 shows the variation 
of capacitance of the typical capacitor cell 
C as a function of voltammetric cycles 
between -800 mV and 800 mV at scan rate 
100 mV sec* 1 . After a slight decrease for 
initial few cycles, the values of capacitance 
become almost constant up to 1000 cycles 
even more. The initial decrease in 
capacitance is possibly attributed to the 
charge consumption in some possible 
irreversible faradic reaction(s) associated 
with loosely bound surface groups at the 
electrode-electrolyte interfaces. 



Fig. 10- Capacitance of the capacitor cell C : ACF | 
PVdF-HFP (20 wt %) - EC:PC (1:1 V/V)-1.0 
M TEABF 4 | ACF as a function of 
voltammetric cycles at the scan rate 100 mV 
see' 1 (Cell area : 2.0 x 1.8 cm 2 ). 
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Polypyrrole and poly(3-methyl thiophene) based 
solid-state redox supercapacitors 

Polypyrrole (pPy) and poly (3 -methyl 
thiophene) (pMeT) based type I and type-II 
solid state redox supercapacitors have been 
constructed by us using polymeric 
electrolytes, PVA-H3PO4 blend and PEO- 
NaC104 complex plasticised with PEG and 
gel electrolyte, PMMA-EC-PC-NaClC^. A 
comparative study on the capacitors with 
symmetical type-I and asymmetrical type-II 
cell configurations have been carried out. 
Results based on galvanostatic charge- 
discharge method have been discussed. 

The polypyrrole and poly (3 -methyl 
thiophene) were electrochemically deposi- 
ted on indium tin oxide (ITO, Balzers, 
sheet resistance ~80 ohm cm 2 ) coated glass 
substrates in a one compartment three 
electrode cell, following the procedure 
described by us elsewhere. 24 Different 
model capacitor cells have been constructed 
using p-doped polypyrrole and poly (3- 
methyl thiophene) with different polymer / 
gel electrolytes. Their configurations are 
reported in Table 5. 

Charge-Discharge studies 

The galvanostatic charge-discharge 
experiments have been carried out for all 
the capacitor cells. The typical charge- 
discharge characteristics of the PMMA-gel 
electrolyte based capacitor cells are 
illustrated in Fig. 1 1 (a-d). The discharge 
capacitance values were calculated from the 
linear portion of the discharge 
characteristics of different cells [Figll (a- 
d]. The values of internal resistance 
R in , (evaluated from initial sudden voltage 
drop in discharge curves)) discharge 
capacitance 'C' (evaluated from eq. 9) and 
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coulombic efficiency 'rf (evaluated from The pPy based symmetrical type-I cells 

eq. 10) of different capacitor cells are listed were charged up to 1.0 volt, which was 

in Table 5. determined by the potential range over 




Fig. 1 1- Galvanostatic charge-discharge characteristics of different capacitor cells; (a) pPy | PMMA-EC/PC-NaC'104 | pPy; 
and (b) pMeT | PMMA-EC/PC-NaC104 | pMeT, (c) (-)pPy | PMMA-EC/PC-NaG104 | pMeT (+) ; and (d) (+) pPy | 
PMMA-EC/PC-NaC104 | pMeT (+). (+)/(-) signs show the polarity of the applied current for charging the 
capacitor cells (charging current = 100 DA and cell area- 4.0 cm-2). 

Table 5-Charge-discharge characteristics of different polypyrrole and poly (3-methyl thiophene) based redox capacitor 


cells. 


Cell configuration 

Type 

Working 

voltage 

(V) 

Internal 

Resistance 

Ri (Kohm cm 2 ) 

Capacitance 
mF cm * 2 

Coul. 

Eff. 

n (%) 

pPy | PVA-H 3 PO 4 1 pPy 

I 

- 1.0 

0.9 

10.0 

90 

pMeT | PVA-H 3 PO 4 1 pMeT 

I 

~0.8 

40.0 

2.5 

82 

pPyj PEO~NaCI0 4 -PEG | pPy 

I 

- 1.0 

36.0 

8.0 

84 

(-) pPy | PEO-NaC10 4 -PEG | pMeT (+) 

II 

- 1.2 

- 

10.0 

- 

(+) pPy | PEO-NaCI0 4 -PEG | pMeT (-) 

II 

<1.0 

8.0 

2.2 

- 

pPy | PMMA-EC/PC-NaC10 4 1 pPy 

I 

-1.0 

1.5 

15.0 

70 

PMeT | PMMA-EC/PC-NaC10 4 | pMeT 

I 

-0.5 

- 

10.0 

78 

(-)pPy | PMMA-EC/PC-NaC10 4 1 pMeT (+) 

II 

- 1.1 

1.4 

12.0 

64 

(+) pPy | P MM A-EC/PC-N aC 10 4 1 pMeT (-) 

II 

-0.6 

1.0 

8.0 

- 
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which the pPy electrode can be doped and 
dedoped. These cells show almost 
symmetrical charge-discharge characteris- 
tics [Fig. 10a] which are indicative of 
identical and symmetrical electrode- 
electrolyte interfaces. The pMeT based 
type-I cell with PVA-H3PO4 blend also 
shows the similar behaviour. While 
charging the pMeT electrode and gel 
electrolyte based type-1 capacitor cell a 
Steep rise in current after about 0.5 V was 
observed [Fig. 10b]. This indicates the 
limitation of the working voltage of the cell 
upto 0.5 V. 

The relatively higher voltage (>1.0 volt) 
is expected in the type-II asymmetric 
pPy/pMeT cells due to different voltage 
range of doping-dedoping of pPy and pMeT 
electrodes. 15 The asymmetric type-II cells 
show peculiar characteristics dependent 
upon polarity of applied charging current. 
The cell configuration (+)pPy | PEO- 
NaC104 -PEG | pMeT(-) shows steep rise in 
voltage after - 1.0 volt which limits the 
range of working voltage <1.0 volt. Also, 
the linear capacitive region has been 
observed for a shorter range of voltage. In 
contrast, the cell with reversed polarity 
functions for slightly higher voltage (~1.2 
volt) and linear capacitive region has been 
observed for relatively larger voltage range. 
This asymmetric polarity dependent feature 
and the capacitance values have distinctly 
been observed in the gel electrolyte based 
capacitors [Fig. 10c&d]. It appears from the 
observed feature that one of the interfaces 
in the asymmetric configuration of the 
capacitor cells becomes dominant over the 
other in terms of charge retention. It further 
suggests that the redox reactions at the two 
interfaces are not truly reversible. In order 
to confirm this aspect, some more 
electrochemical studies would be required. 
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Redox Supercapacitors based on MnC>2- 
Polypyrrole Composite Electrodes 

Redox Supercapacitors using electro- 
chemically synthesized MnC>2-polypyrrole 
composite electrodes have recently been 
fabricated by us with different electrolytes, 
namely polymer electrolyte film, PVA- 
H3PO4 aqueous blend, aprotic liquid 
electrolyte, LiClC>4-PC and polymeric gel 
electrolyte, PMMA)-EC-PC-NaC10 4 . 

MnC>2 is well known to be an important 
material used in primary dry cells as de- 
polariser and in rechargeable lithium 
batteries as active cathode. 32 ' 34 This aspect 
of MnC>2 was considered to be exploited as 
a possible promising electrode material for 
supercapacitors. However, since MnC>2 
exhibits poor conductivity, a conducting 
network is essentially needed to use it as an 
electro-active electrode material. Generally 
activated carbon or acetylene black powder 
is dispersed in MnC>2 to make it conducting. 
Recently, S. Kuwabata et al 35 reported a 
composite MnCVpolypyrrole (MnCh-pPy) 
for its application in lithium secondary 
batteries as cathode material. The MnC>2- 
pPy composite works as conducting 
network of electro-active material in thin 
film form. Efforts have been made to utilize 
this conducting network as electro-active 
electrodes for the fabrication of solid state 
redox supercapacitors. 

Secondly, polypyrrole itself is a poten- 
tial conducting polymer electrode for 
secondary batteries and Supercapacitors. 1 ' 3, 
| s- | 9, 22 Yh e utilisation of the polypyrrole 
can be enhanced by preparing the highly 
porous film which may permit the easy 
insertion-deinsertion of ions at the elec- 
trode-electrolyte interfaces. Consequen- 
tly, the enhancement in the energy/power 
density is possible with the porous 
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polypyrrole films. Many authors reported 
different ways to enhance the porosity of 
polypyrrole/ 6 The Mn0 2 -pPy system also 
acts as a porous composite film, which may 
have an added advantage for its possible 
application as electrodes in redox super- 
capacitors. 

In view of the above facts, the Mn0 2 - 
pPy composite electrode material has been 
used to fabricate redox-supercapacitors. A 
comparative study has been carried out on 
the supercapacitors with different electroly- 
tes, solid polymeric electrolyte, PVA- 
H3PO4 blend, aprotic liquid electrolyte 
LiC104-PC and sodium ion conducting gel 
electrolyte, PMMA-EC/PC-NaC10 4 . The 
performance characteristics of the capa- 
citors have been evaluated by galvanostatic 
charge-discharge tests. 

The Mn0 2 -pPy composite electrodes 
were electrochemically deposited on ITO 
coated glass substrates. The electrode- 
position was carried out at constant current 
in single compartment three electrodes cell. 
The procedure of electrosynthesis has been 
described in detail by us elsewhere. 25 
Following three model capacitor cells have 
been constructed using Mn0 2 -pPy com- 
posite electrodes and three different 
electrolytes : 

Cell 1 : Mn0 2 -pPy I PVA-H 3 P0 4 | Mn0 2 - 
pPy 

Cell 2 : Mn0 2 -pPy | UC10 4 -PC | Mn0 2 - 
pPy 


Cell 3 : Mn0 2 -pPy | PMMA-EC/PC- 
NaC10 4 | Mn0 2 - pPy 
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Charge-Discharge studies 

The galvanostatic charge - discharge 
characteristics of the capacitor cells 1 -3 are 
shown in Fig. 12. All the cells were 
charged up to 1.0 volt, which is fixed 
within the voltage range of electrochemical 
stability of the electrolytes. The 
symmetrical polypyrrole based capacitors 
are usually charged up to a 1.0 volt limit 
which is determined by the potential range 
over which the polypyrrole based 
electrodes can be doped and dedoped. 

The typical charge-discharge character- 
istics (i.e. voltage vs. time plot) of the cell 
Mn0 2 -PPy | PVA-H3PO4 | Mn0 2 -pPy 
shows a non-linear behaviour [Fig. 12(a)], 
This indicates a non-capacitive behaviour 
of the cell. The charge-discharge 
characteristics has been compared with the 
similar cell with pure polypyrrole 
electrodes, as shown in Fig. 12(b). This 
shows a typical linear discharge 
characteristics which is an indication of 
capacitive behaviour of the cell. The non- 
linear behaviour of the charge-discharge 
characteristics confirms the role of Mn0 2 
present in the composite network of Mn0 2 - 
pPy system. This is possibly owing to the 
following reversible chemical reaction of 
Mn0 2 with water based proton conducting 
electrolyte, PVA-H3PO4 blend. 37 

Mn0 2 + H 2 0 + e ► MnOOH + OH' 

In order to confirm this aspect, a 
capacitor cell has been constructed with 
aprotic liquid electrolyte LiC10 4 -PC to 
avoid the above aqueous based reaction. 
The discharge characteristics of the cell 
(Fig. 12c) is almost linear which indicates 
its capacitive behaviour. This is possibly 
owing to the easy insertion and de-insertion 
of Li + ion through the Mn0 2 matrix along 
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with the conventional Faradaic charge- 
discharge reaction with polypyrrole, 
leading to the pseudo-capacitance. This 
comparison indicates that the MnC^-pPy 
composite is a suitable electrode system 
only for aprotic (non-aqueous) electrolytes. 


shows the typical variation of discharge 
capacitance with respect to the varying 
current density for the cell Mn02-pPy | Gel 
Electrolyte | Mn02-pPy. The stable values 
of capacitance (10.0-18.0 mF cm' 2 ) for 
different discharge current densities have 
been observed. The internal resistance of 






Fig. 12- Galvanosiatic charge-discharge characteristics of different capacitor cells; (a) Mn0 2 -pPy | PVA-H 3 P0 4 1 
MnO?-pPy, (b) pPy | PVA-H 3 P0 4 1 pPy, (c) Mn0 2 -pPy | LiC10 4 -PC | Mn0 2 -pPy, and (d) MnO r pPy | 

PMMA-EC/PC-NaC104 | MnO r pPy. 


A solid state redox capacitor cell based 
on Na + ion conducting gel electrolyte 
(PMMA-EC/PC-NaC10 4 ) has also been 
fabricated which shows a capacitive 
charge-discharge behaviour, as shown in 
Fig. 12(d). The values of discharge 
capacitance (evaluated from eqn. 9), 
internal resistance Rj (evaluated from initial 
sudden voltage drop in discharge curves) 
and the coulombic efficiency, rj, (evaluated 
from eqn. 10) are listed in Table 6. Fig. 13 
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the cell (observed to be -1500 ohm cm 2 ) 
has been compared to that of the cell with 
pure polypyrrole electrodes and almost 
same value of internal resistance has been 
observed. This indicates that the 
incorporation of MnC >2 in polypyrrole 
exhibit no substantial effect on the internal 
resistance of the capacitor cell. The 
polypyrrole plays the effective role to form 
conducting network in the MnCVpPy 
composite system. Finally, it has been 
concluded that the MnC> 2 -pPy composite 
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Table 6- Charge-dischargecharacteristics of different Mn0 2 - 

-polypyrrole based redox capacitor cells. 


Cell configuration 

Internal 
Resistance Rj 
(ohm cm 2 ) 

Discharge Capacitance 

C (mF cm' 2 ) 

Coulomb 

ic 

Efficienc 

y n (%) 

Mn0 2 -pPy | PVA-H 3 PO 4 1 Mn0 2 -pPy 

5xl0 4 

Non-capacitive 

- 

MnO r pPy 1 LiC10 4 -PC 1 Mn0 2 -pPy 

1.5xl0 3 

- 10.0 

44 

Mn0 2 -pPy | PMMA-EC/PC-NaC10 4 1 Mn0 2 -pPy 

1.5xl0 3 

10.0-18.0 

66 


system has been proved to be a potential 
electrode material for fabrication of redox 
supercapacitors using aprotic (non-aquous) 
liquid or polymer/gel electrolytes. 



Fig. 13- Variation of discharge capacitance of the cell 3: 
Mn0 2 -pPy I PMMA-EC/PC-NaC10 4 | Mn0 2 - 
pPy as a function of current density. 


Applications 

There has been increasing global 
interest in the development of 
‘ Supercapacitors ’ due to their various 
possible technological applications as 
power source, mentioned below: 

• Portable electronic equipment 

• Computer memory back-up 

• Medical equipment 
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• Power supplies 

• Load levelling 

• Electrical Vehicles 

In low energy density applications such 
as electronic/medical devices, supercapa- 
citors may be used as an alternative power 
source, whereas it is a complementary 
power source to the batteries in electrical 
vehicle like high energy density applica- 
tions 1 ' 3 . One of the most promising 
application of Supercapacitors is expected 
in the hybrid power design for electrical 
vehicles where they would be used as 
complementary power system to the main 
rechargeable batteries. The battery system 
would provide energy for continuous 
motion where the Supercapacitors (connec- 
ted in parallel to the batteries) ensure peak 
power delivery for acceleration. 
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The National Academy of Sciences, 
India, founded in 1930 with Prof Meghnad 
Saha, one of the doyens of Indian Science 
in the early twentieth century, as its first 
President, is the oldest science academy in 
the country. The main objective of the 
Academy is to provide a national forum for 
the publication of research work carried out 
by Indian scientists and to provide 
opportunities for exchange of views among 
them. The Academy, starting with 19 
fellows and 57 ordinary members grown to 
its present size of nearly 1Q12 fellows and 
1350 members, is taking a leading role in 
fostering research in frontier areas of 
sciences by arranging topical symposiums 
during its Annual Sessions. The Academy 
is also trying to mould national consensus 
on important issues of both scientific and 
societal interest by holding discussion 
sessions and symposiums on issues related 
to science and society that have an 
important bearing in the national context. 
An important activity of the academy is to 
provide opportunity to young scientists to 
participate in the Annual Sessions and 
present their research work. This allows the 
younger generation to interact with peers in 
their own field and benefit from the 
discussions with their peers. The peers, in 


turn, have the opportunity of getting 
exposed to some radically new ideas, which 
can only come from the younger scientists, 
who are not burdened with preconceived 
ideas. It is indeed a laudable activity that 
will go a long way in fostering appropriate 
research culture among the upcoming 
young scientists. 

The lack of effort to impart science 
education in a manner that will generate 
curiosity and genuine interest and nourish 
creativity and innovativeness among the 
younger generation in the country has been 
a matter of grave concern. This led the 
Academy to start a programme of Science 
Communication during the last decade 
involving organisation of various activities 
for school students and refresher courses 
for college teachers. This programme has 
now been widened in its scope as a result of 
the Vision Paper prepared during the 
current year, to include intensive 
discussions on topical issues of Science and 
Technology that have direct societal 
relevance. 

The Current Session 

The Annual Session of the Academy is 
being held at Ahmedabad after nearly a 


* Based on the inaugural address at the 73 rd Annual Session of The National Academy of Sciences, India, delivered 
on Oct 10, 2003 at Physical Research Laboratory (PRL), Ahmedabad. 
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quarter century and for the first time at 
PRL. The aim of the Academy to foster 
relationship with as large an academic 
community as possible during such a meet 
is also well served by the fact that the 
Physical Research Laboratory and the 
Gujarat University are co-hosting this 
Session with the University holding the 
Biological Sciences Session in its premises. 
The interest shown by the academic 
community in participating in the 
deliberations of the scientific sessions is 
evident from the presentation of nearly 200 
papers. The gesture of the Academy and the 
host institute to provide support to a large 
number of young students and research 
scholars from Gujarat and the neighbouring 
States to participate in the Annual Session 
is indeed very laudable. 

The Physical Research Laboratory has a 
very strong tradition in the field of 
astronomy and astrophysics. PRL is unique 
in that, under one umbrella, it has been 
studying physical phenomena related to the 
earth's interior core, upper atmosphere, 
ionosphere, magnetosphere, interplanetary 
space, cosmic ray propagation, stellar 
astronomy and cosmology. It was in 1936 
Prof Saha published a paper titled "On a 
stratospheric astrophysical laboratory" 
bringing out the advantage of carrying out 
observations in space. PRL with its 
pioneering work in the fields of cosmic 
rays, X-ray astronomy and radio-astronomy 
realised the importance of space based 
observations, thanks to Dr Vikram 
Sarabhai's Vision in the fifties and started 
building up its space capability. PRL, under 
the leadership of Dr Sarabhai, thus became 
the cradle of the Indian Space Programme. 
Starting with sounding rocket experiments 
from Thumba for investigating the upper 
atmosphere, equatorial electrojet, and X-ray 
astronomy phenomena, PRL was respon- 
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sible for initiating steps to build a strong 
space science and technology group. At 
present PRL has active research pro- 
grammes in the fields of infra-red and 
optical astronomy carried out from the 
Observatories at Mt. Abu and Udaipur 
respectively. A low energy X-ray payload, 
designed and developed at PRL, that was 
launched into space early this year by 
GSLV, is currently collecting valuable data 
on solar X-rays which have already yielded 
interesting results. These can shed new 
light on our understanding of the solar flare 
phenomena over a wide scale, from major 
flares to nano-flares. Laboratory studies 
that have important bearing on solar system 
astrophysics are also being pursued 
actively. Plans are afoot to enter into the 
new field of sub-millimeter astronomy and 
enhance the capabilities of the solar 
observatory at Udaipur. PRL is also taking 
a lead role in ISRO's ambitious future 
programme on planetary exploration with 
the proposed Lunar Polar Orbiter mission 
in 2008 as its forerunner. Apart from the 
experimental initiatives, PRL hosts a very 
strong theoretical group in Astro-particle 
physics and cosmology that has made very 
significant contributions, particularly in the 
field of neutrino astronomy, a very active 
area of current research in the global 
context. It is therefore very appropriate that 
a national symposium on "Astronomy in 
the New Millennium" is being held as a 
part of this Annual Session of the 
Academy. 

Observational astronomy has made 
rapid strides during the last decade both in 
the global and national context. Space- 
based platforms, such as the Hubble 
Telescope, the "Chandra" Observatory, and 
the recently launched Space Infrared 
Telescope Facility (SIRTF) offer new vistas 
in observational astronomy, which have led 
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to significant refinement in our understand- 
ing of the early evolution of the universe 
and formation and evolution of various 
stellar objects. Significant improvements in 
performances of ground-based astronomical 
Observatories as well as establishment of 
new facilities across the globe have 
complimented the space-based observations 
and consolidated our understanding of 
astrophysical processes both at small and 
large spatial and temporal scales. 

India has a long-standing tradition in 
observational astronomy starting with the 
study of solar phenomena from the 
Kodaikanal Observatory in the early 
twentieth century to the recent develop- 
ments such as the establishment of the 
Giant Meter Wave Radio Telescope 
(GMRT) near Pune, the high altitude 
observatory at Hanley for observations in 
the optical and infrared regions and the 
gamma-ray Observatory at Mt. Abu. There 
are ongoing plans for several new 
Observatories as well as for upgradation 
and expansion of existing facilities to 
venture into new research areas in the field 
of solar, stellar and extragalactic 
astronomy. New initiatives for space-bome 
astronomical observations using indigenous 
satellite launch capabilities have been 
successfully implemented and work has 
been initiated for a dedicated astronomical 
satellite, the Astrosat. An ambitious 
programme for the exploration of our own 
solar system is also being launched. 

Research in the field of Astronomy and 
Astrophysics in India is thus poised for a 
major advance and the proposed 
symposium on "Astronomy in the new 
millennium" is extremely timely. A broad 
spectrum of subjects (cosmology and early 
universe, evolution of stars and galaxies, 
high energy stellar phenomena, neutrino 
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astronomy, X-ray and gamma ray 
astronomy, Sun and solar system studies) 
are being covered during the symposium by 
leading experts from various Institutions 
and Observatories. It is hoped that the 
Symposium will review critically the global 
vs local status in this field of research that 
will lead to a better definition and focus of 
the new programmes and at the same time 
provide an overall perspective of this 
exciting field of research to the broad 
spectrum of people from various branches 
of sciences constituting the fellowship and 
membership of the Academy. 

Role of Academy 

I thought I should touch upon some 
issues, which in my opinion are critical to 
the future functioning of this Academy as 
well as other similar Academies. I believe 
that the academy should be a think-tank 
and a pressure group. It should initiate and 
encourage open discussions on topical 
issues of science and technology and on 
issues of societal importance. A few 
suggestions are made here for the 
consideration of the Academy: 

1. Induction of young scientists - Only 12 
out of 1012 Fellows of the academy are 
below the age of 45. The proportion of 
younger scientists in other academies is 
still worse. This does not auger well for 
the future growth of science in India. A 
concerted effort by the fellows is 
required to identify and induct bright 
young scientists in India as fellows of 
the academy. 

2. Very little is done in promoting science 
and technology education in the 
country, which is in a very bad shape. 

• Revitalisation of education by providing 
training to teachers, promoting quality 
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research and encouraging younger 
generation to share the excitement of 
science are some areas in which the 
Academy can play a positive role. 
Massive use of distance education is 
necessary to accomplish these aims. 

We have a catch 22 situation in the 
publication of journals. Good papers are 
not published in our journals, which 
results in poor citation. Poor citation 
does not encourage scientists to publish 
good papers. It may be worthwhile to 
have journals carry only expert review 
articles on various subjects along with 
short papers/notes. 

Ethics of Science 

Since science and technology deal with 
people's requirement for better living, it 
must necessarily be based on ethical 
principles. Historically two lines of thought 
emerged, the first derived from Plato and 
Descartes and the second by Aristotle. The 
Plato approach (principle based) based on 
universal principles derived from pure 
reason without regard to the situation, 
culture or consequences, was later modified 
by Barker and Mill to make it consequence 
oriented. This resulted in the moral choice, 
if we may call it so, being bestowed on 
utilitarian yield - the greatest good for the 
greatest numbers. Aristotle's ideas, 
(practical wisdom), which were forgotten 
for a considerable time, have been 
resurrected by contemporary philosophers 
such as Dreyfus to take into consideration 
the variable of specific situations, practices 
and norms. 

Application of ethics, however, have 
got very complicated in the modern age 
because of the range and influence of 
scientific actions, which can have local 
(building collapse, shuttle explosion), large 
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area (environmental) or even global effect, 
(green house emission) and manifest very 
quickly (over the internet), very slowly 
(radio-active disposal), temporarily persis- 
tent (organic pollutants) or become even 
irreversible (genetically modified organism 
and cloning). It is therefore all the more 
necessary that scientists integrate societal 
and technical parameters in their thinking 
not only for the development but also for 
the survival of humankind. 

In order to promote the above objective, 
the Academy must become the conscience 
of scientific community. Penchant for going 
to the press without going through peer 
reviews has caused a lot of damage. Cold 
fusion, plagiarizing, sensationalisation 
(recent controversies on pesticide/chemical 
contamination of cold drinks) are some 
examples of the above trend. Pseudo self- 
styled environmental scientists making 
statements, starting agitations merely for 
publicity (environment, dams, energy, 
pollution, etc.), without caring for the 
consequences can only cause great damage 
to societal development. On the other hand 
the habit of scientists and engineers to 
promote large programs, without studying 
their consequences, is equally regrettable. 
Talking on large dams like Hirakund dam 
and Damodar valley project, Prof Saha, in 
1954, chided the scientific community of 
big dam enthusiasts with these words "we 
shall make the highest dams in the world, 
may be a very good slogan, but is very bad 
from the scientific point of view". Instead 
of taking our decisions on objective cost- 
benefit analysis, we often seem to give our 
advise palatable to the whims of the self 
seeking politicians. 

Academy should consider setting up of 
an ethics committee of eminent scientists to 
look into matters, which relate to misuse, 
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sensationalisation and even general conduct 
of scientists. Larger issues such as, ethics 
of xeno-transplantation, cloning, energy, 
space mines, destablisation of food chain 
by indiscriminate activities etc., need to be 
discussed objectively and dispassionately. 
Multi-national companies with their 
monopoly of terminator seeds can control 
world's food supply. Infection of fields of 
farmers, who have not used this 
technology, by pollen from crops carrying 
the terminator gene is a real problem. 
While the crops in the particular season 
will not be affected, some of the seeds may 
become sterile when they are sowed in the 
next season, affecting the entire food chain. 
Use of space mines, which is being pushed 
in the name of security, can itself 
accidentally wipe out the space assets of 
many nations. Unless Academies take a 
proactive role in matters which affect the 
society as a whole, I believe they will 
remain as elite clubs at best, far removed 


from reality. Intensive subject discussions 
can and are best held during purely 
scientific meetings/workshops in particular 
identified areas. Larger issues can only be 
discussed and dealt with by the academies. 

I wish to conclude with the statement 
made by the great scientist Peter Kapitza on 
Soviet Academy, which applies equally to 
our Academy, and which was often quoted 
by Prof M N Saha: "The academy of 
science is the general headquarters of 
science and should be the idealogical guide 
of the science of the whole nation. Every 
individual institute and member of the 
academy must follow the same policy, 
strive to become the guiding influence in 
their particular sphere of science and try to 
achieve the maximum development". I 
hope this academy will take greater interest 
in harmonizing scientific progress with 
societal aspirations. 
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Abstract 

Urease (3, 5. 1.5) from jack bean meal was 
immobilized by embedding in dry films of 
reinforced bees wax and by covalent linkage using 
ethylene diamine-glutaraldehyde arm onto seed 
coats (pellicles) of garden cress {Lepidium sativum) 
in aqueous suspension. The Activity of free and 
immobilized enzyme as a function of pH, 
temperature, storage stability, kinetic parameters 
and periodic use were compared. Enzyme bound 
under mild conditions retained good activity and 
storage stability during repeated use as compared 
to the free urease. Kinetic parameters viz K m and 
V max are in accordance with Michaelis-Menten 
rate equation 

(Keywords: urease/seed coats/beeswax/immobili- 
zed enzyme) 


Natural supports offer ideal means for 
immobilizing enzymes. A number of 
enzymes have been immobilized using 
temperate conditions on various natural 
substances such as seed coats, cellulose 
etc. ’ . These are extremely cheap, eco- 
friendly and ready to use as compared to 
synthetic supports. Amongst supports used 
for immobilization, biological macromole- 
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cules like cellulose and its derivatives 
through physical adsorption or chemical 
linkages have been studied most 
extensively 3,4 . Urease has been immobili- 
zed on many supports like calcium alginate 
paraffin wax, lac and some natural seed 
coats 5,6 . Immobilization of enzymes on 
beeswax, a natural animal product is yet not 
reported. Bees wax is a mixture of myricyl 
palmitate [CH3-(CH2)i4-COOC3oH6i] and 
n-hexacosanyl palmitate [CH 3 -(CH 2 )i 4 - 
COO-C26H53]. This substance finds enor- 
mous applications in food, pharmaceutical, 
cosmetics, microencapsulation, polishes 
and leather industries. 

When enzymes are bound under mild 
conditions on polysaccharides activated by 
various means, mechanical instability is a 
major problem. Another difficulty that 
arises is due to their particle size, which 
makes them readily compact resulting in 
high pressure drops in packed bed column 
reactors. To overcome this problem, the use 
of pellicular catalysts, which consist of 
enzyme immobilized only on the outer 
porous shells of carrier particles, is highly 
recommended 7 . Because of its more 
exposed outer structure within an inert hard 
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core, Ocimum basilicum seeds have served 
as efficient natural pellicular support for 
immobilizing enzymes 1 . 

Lepidium. sativum seeds are used as 
diuretic, tonic, aphrodisiac, laxative, 
rubifacient in poultices for hurts and 
sprains. Seeds also yield a semidrying oil, 
used for burning and soap making. The 
seed mucilage called 'cress seed mucilage’ 
is used as a substitute for tragacanth and 
arabic gum. It allays irritation of the 
intestines in diarrhoea and dysentry. The 
present paper reports the use of alternative 
Lepidium sativum seeds, which are 
mucilaginous, nontoxic [edible] and offer a 
cheaper and natural polysaccharide support 
for enzyme immobilization. Also, same 
enzyme is immobilized on bees wax, a non- 
toxic biocompatible substance. 

The enzyme urease from jackbean meal 
(Loba Chem. Pvt. Ltd. India.), -its substrate 
urea (Loba ) and glutaraldehyde (Loba) 
were used and 1,4- butenediol diglycidyl 
ether was a Sisco product. Nessler's 
reagent, Tris buffer and other chemicals 
used were of Analar grade. Muslin cloth 
(Tata fabric pure cotton) was used. 
Lepidium sativum seeds were procured 
from local market and beeswax was 
extracted from honeycomb as described 
below. 

Dry honey combs of Apis dorsata were 
chopped into smaller pieces and boiled in 
hot water in an aluminium vessel for 10 - 
15 minutes till all the cells of honey combs 
separated from each other and yellow 
coloured melted wax started to float on the 
surface. This was then strained through a 
strainer made of fine mosquito net cloth. 
The filtrate contained water plus wax 
floating on the surface, which cooled to 
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solid. The wax layer from the top thus 
collected was again boiled in distilled water 
and wax, free from insoluble material was 
collected. The wax thus obtained was 
further- purified by washing in methanol 
first and then in acetone and was used for 
further experiments. 

(1) Preparation of immobilized urease on 
bees wax films: 

To 4g of melted beeswax in 
thermostated waterbath. lg of enzyme 
powder was added and kept on stirring 
continuously at slow speed. Muslin cloth 
was previously washed in distilled water 
and sundried. Of this 1 sq.cm pieces were 
cut (wt. 5mg.). Each of these pieces was 
dipped in wax enzyme mixture which was 
continuously agitated and kept in water 
bath thermostated at 60°C. After 3-5 
seconds these pieces were removed by 
using pointed B.B. forceps, slightly shaken 
to remove dripping wax and allowed to dry. 
In this way 1 cm 2 reinforced film of 
beeswax with enzyme was made. Films 
were weighed to find out amount of 
enzyme immobilized. 

(2) Preparation of immobilized enzyme on 
seed coats: 

Dry seeds of Lepidium sativum were 
swollen by soaking in distilled water at 
room temperature for 20-24 hours. These 
were then epoxy activated by method of 
Sundbfrg and Porath 8 after slight 
modification. To 20 ml of swollen seeds 
1 0ml of 1 -4 butenediol diglycidyl ether and 
10 ml of 0.6 N NaOH was added and kept 
on stirring for 4 hours and then they were 
left undisturbed for another 20 hours. After 
this- reaction was stopped by washing the 
seeds with ample water. These epoxy- 

SCIENCE LETTERS, VOL. 27, NO. 1 & 2, 2004 



SHORT RESEARCH COMMUNICATIONS 


activated seeds were then treated with 10 
ml of 0.1M ethylene diamine (pH 11.0), at 
room temperature for 24 hours with 
occasional shaking. To this reaction 
mixture, 10 ml of 2% glutaraldehyde was 
added and kept for another 24 hours with 
initial shaking only. These were then 
washed with excess of water. 

A solution of urease (5mg/ml ) was 
added to gluteraldehyde activated seeds and 
kept in cold with intermittent shaking for 
24 hours. The reaction was terminated by 
washing the seeds with distilled water and 
then with phosphate buffer pH 7.0. Seeds 
with immobilized enzyme were stored in 
phosphate buffer solution in glass 
stoppered bottles at 4 °C. 

Free enzyme in the washings was 
measured colorimetrically using Lowry's 
method to obtain the amount of bound 
enzyme. The amount of enzyme 
immobilized on the support was calculated 
later. 

Enzyme assays 

(1) Assay of native urease : The 
standard solution was prepared by 
dissolving lOmg urease in 10ml phosphate 
butler (pH 7). 1ml of 3% urea solution and 
1 mg/ml of enzyme was taken and 
incubated for 15 minutes at pH 7 and 
temperature 55°C in a water bath, The 
reaction was stopped by placing the 
reaction mixture in ice-bath and 0.66N 
H 2 SO 4 was added. 1ml of Nessler's reagent 
was then added to develop colour and 
optical density of the solution, was 
measured using ESICO colorimeter at 500 
nm. 

(2) Assay of urease immobilized on 
beeswax film : One film of beeswax was 


taken in 1 ml buffer solution and reacted 
with 1 ml 3 % urea solution for 15 minutes 
in water bath. The same procedure followed 
the reaction and measurements were similar 
to the native enzyme except that the films 
were taken out from the reaction mixture 
just before stopping the reaction. 

(3) Assay of enzyme immobilized on 
seed coats : To 1ml of seeds containing 
immobilized enzyme 1 ml of phosphate 
buffer and 1 ml 3% urea solution was 
reacted for 1 5 minutes. The same procedure 
followed the reaction and measurements 
were similar to that carried out for native 
enzyme except that seeds were removed 
before chilling the solution and adding 
H2SO4. 

Determination of enzyme activity 

The catalytic activity of native and 
immobilized enzymes was measured 
colorimetrically at 500nm. After reaction 
with the substrate urea, the amount of 
ammonia liberated was measured using 
Nessler's reagent method • as described 
before. 

Determination of kinetic parameters 

For determination of kinetic parameters 
substrate concentration was kept on 
changing keeping enzyme concentration 
fixed for both native as well as 
immobilized enzyme system at constant 
temperature 40°C and pH=7.0. 
Measurements were made as described 
before using Nessler's reagent method. 
Absorbance was read at 500 nm. Velocity 
was measured in moles of NH 3 produced / 
min / mg enzyme. K m and V max were 
determined using Lineweaver-Burke plot. 
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Results and Discussion : Due to mild 
and precise conditions of physical and 
chemical treatment during the immobilizing 
procedure the structural integrity and 
activity of urease was preserved. The 
entrapment of urease in the beeswax was 
such that it left some part of enzyme 
exposed to carry out catalytic function and 
in such a condition it retained 40% activity. 
Urease is a thermostable enzyme that is 
denatured only at 80°C and therefore it 
could easily be immobilized on beeswax at 
60°C. The reinforced films of bees wax 
were durable. 

L. sativum seeds like other mucilagi- 
nous seeds swell upon wetting. Seed 
mucilage known as cress seed mucilage is 
used as a substitute of tragacanth gum and 
is used to allay irritation of intestines in 
dysentery and diarrhoea 9 . This pectinous 
outer layer substance is termed as amyloid 
and is a galactosan [polysaccharide] having 
empirical formulae [CeHio Os],, 10 . It has a 
large capacity for hydration. One gram of 
dry seeds on soaking result in about 15 ml 
of swollen seeds. The mucilage consists of 
thread like microfibrillar structure with a 
large surface area as observed under a 
microscope. The seed consisting of a hard 
core with pectinous porous outer layer can 
serve as a cheap and ready' to use natural 
polysaccharide support for immobilizing 
enzymes. The polysaccharide layer can be 
modified by epoxy activation for successful 
binding of urease using covalent attachment 
via. ethylene diamine-glutaraldehyde arm. 
Enzyme thus bound (2mg/ml seeds) 
retained 63% activity. 

No leakage of enzyme was observed 
after continuous washing the preparation 
with 0.2 M phosphate buffer or distilled 
water. The column of urease bound seeds 
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maintained high flow rate [10ml column/30 
minutes for 10 ml buffer]. Studies were 
carried out with these preparations and 
compared with native enzyme to find out 
utility of urease in hydrolyzing urea. 

pH activity profile 

The pH dependence of the relative 
activity of the enzyme urease is compared 
with that of the enzyme immobilized on 
various supports. The reaction was carried 
out in 0.2M phosphate butter for pH range 
5.6-7.8. The optimum pH for native and 
immobilized enzymes on seeds were 6.8 
and 7.0, respectively. This shows that 
immobilization shifts the optimum pH 
range to basic side by 0.2 units but in the 
case of enzyme immobilized on bees wax 
film showed optimum pH at 7.2, a shift of 
0.4 units to basic range. This shows urease 
can be used in narrow basic range than 
native urease. (Fig. 1) 



pH 

Fig. 1- Plot of percentage maximum activity vs. pH 
for pH profile of free (•) and immobilized 
urease on beeswax (A) and seeds (■). 

Thermal stability 

The optimum temperature for free as 
well as immobilized enzyme was around 
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55 °C, when thermal stability studies were 
noted at different temperatures and at 
constant pH=7 kept for 15 minutes. Both 
the native as well as immobilized enzyme 
were stable upto 75 °C. However, the 
activity of immobilized enzyme on seeds 
was higher at high temperatures as 
compared to native enzyme (Fig. 2). 
Temperature stability studies for enzyme 
immobilized on bees wax could not be 
carried out as the melting point of bees wax 
was around 58°C. 



Temperature <*C) 

Fig. 2- Plot of percentage maximum activity vs. 
temperature for thermal stability of free (•) 
and immobilized urease on seeds (A). 

Storage stability 

The native enzyme retained same 
activity for about 15 days, which slowly 
decreased and was lost completely after a 
month. The enzyme immobilized on 
beeswax retained 98% activity for 15 days 
and 59% activity even after a month in dry 
condition. Enzyme immobilized on seeds 
retained more than 82% activity after about 
15 days when stored at 4°C (Fig.3). 
Moreover no physical change, structural 
alterations or microbial growth was seen on 
seeds stored at room temperature for longer 
periods. 
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Fig, 3- Plot of percentage maximum activity vs. time 
(days) for storage stability of urease 
immobilized on beeswax (•) and seeds (A). 

Retention of activity after reuse 

Immobilized enzyme on beeswax film 
could be used repeatedly upto 13 batches 
under assay condition retaining 90% 
activity while immobilized enzyme on 
seeds retained 60% activity upto 5-6 
batches under assay conditions (Fig. 4). 



Mg 4- Plot of percentage maximum activity vs. 
number of times used for reusability of urease 
immobilized on beeswax (•) and seeds (A). 

Mechanical stability 

Beeswax films were mechanically 
stable when centrifuged or kept on 
magnetic stirrer at full speeds. Seeds were 
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also stable and resistant to destruction when 
stirred continuously for 10-12 hours on 
magnetic stirrer. When seeds were 
centrifuged initially they adhered to each 
other forming a pellet but could be 
separated easily when resuspended in 
aqueous medium. 

Kinetic studies 

The kinetic parameters of K m and V max 
were evaluated using Lineweaver and 
Burke plot (Fig. 5) and are tabulated 
(Table- 1). The kinetic parameters are in 
accordance with Michaelis - Menten 
equation. Immobilization of enzyme leads 
to a decrease in both K m and V max values. 
This shows that active sites of enzyme are 
masked due to immobilization. Further that 
in case of immobilization on seeds K m 
value is further lowered This can be 
attributed to the use of chemicals viz. 
ethylene diamine-glutaraldehyde for 
covalent linkage and therefore more active 
sites are blocked. 



l/S (mol. I -')’ 1 

Fig. 5- Double reciprocal plot of velocity V (moles of 
NH 3 produced/min/mg enzyme.) vs. 
concentration of urea for free (•) and 
immobilized urease on beeswax ( A ) and seeds 
(«). 
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Table 1- Kinetic parameters of free and immobilized 
urease on bees wax and seeds 



K m (moles 

v 

v max 


x lO' 1 ) 

moles/min/ 

mg.enzyme 

Native 

0.250 

1.54 x 10 -4 

Immobilized 
enzyme on beeswax 

0.200 

l.llxlO" 4 

Immobilized 
enzyme on seeds 

0.143 

1.30 x 10 -4 
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Abstract 

Homogeneous solid solutions of cadmium calcium 
hydroxylapatite, containing arsenate, [Can,. 
nCdnfPO^CAsO^mCOH)^ with fixed m = 1 ; n = 

1, 2, 3 have been prepared over the entire 

compositional range by co-precipitation at 37+ .5° 
C in aqueous media by addition of ammonium 
dihydrogen phosphate and ammonium arsenate to 
the solutions of cadmium and calcium nitrate 
made strongly basic with liquor ammonia. The 
infrared spectra and lattice parameters of the 
solid solutions have been measured. 

(Keywords : Hydroxylapatite / solid solutions / co- 
precipitation) 


Introduction 

In recent years, attention has been 
directed towards all the basic calcium 
phosphates for the process of calcification 
in general and for dentistry in particular. 
Hydroxylapatite, represented as [Caio(PC> 4 ) 
6 (OH) 2 ], is now recognized as the most 
important of the inorganic constituents of 
human bone and teeth 1,2 . It belongs to an 
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isomorphous series of compounds known 
as apatite which exists in nature as the 
mineral hydroxylapatite and is similar, if 
not identical, to bone mineral which can be 
prepared in aqueous solution 3,4 . Solid 
solution of cadmium and calcium hydro- 
xylapatites containing arsenate were pre- 
pared separately earlier by heating 
(~1300°C) mixtures containing various 
proportions of cadmium and calcium 
hydroxylapatite with asenate 5 . These 
samples prepared by the solid state reaction 
were, however, found to be non- 
homogeneous. 

An attempt has now been made to 
prepare homogeneous solid solutions of 
cadmium-calcium hydroxylapatite with 
arsenate over the entire compositional 
range by co-precipitation in aqueous 
media 6,7 . 

10-n Ca(N0 3 ) 2 + nCd(N0 3 ) 2 + 6-mNH4H 2 P0 4 + 
m(NH 4 ) 3 As 04 + liquor NH 3 — > 

[Caio- n Cd n (PC> 4 ) 6 . m (As 04 ) m(OH) 2 ] 
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Materials 

Chemicals used for the preparation of 
these samples were of reagent grade. Water 
used in the preparation and in washing was 
boiled to remove CO 2 and then used 
immediately. 

Preparation 

Stoichiometric quantities of ammonium 
dihydrogen phosphate (0.25 M) and 
arsenate (0.25 M) solution A and solution 
B containing (0.25 M) solution of each of 
cadmium and calcium nitrate were prepared 
separately in CO 2 free doubly distilled 
water. Both solutions were initially more 
basic than pH 11 by addition of liquor 
ammonia and this high basicity was 
maintained during the precipitation reaction 
and subsequent digestion. Then a part of 
solution B was put in a flask (2L) fitted 
with two separating funnels and a delivery 
tube. For the preparation of pure CaHA, an 
excess of calcium ion was necessary to 
prevent the uptake of excess phosphate 
accordingly considerable cation excess was 
used throughout. Solutions A and B were 
taken separately in the separating funnels 
and added drop wise to the flask 
simultaneously in a calculated amount. 
Precipitation was carried out by bubbling 
CO 2 free air. The whole was heated under 
reflux for 2 hrs to improve homogeneity 
and crystallinity of the precipitate. Desired 
pH during precipitation was maintained by 
testing the filtrate after separation of the 
precipitate, since any alteration of pH of the 
medium during precipitation leads to the 
formation of cations deficient apatite 8 . The 
precipitate was allowed to settle overnight 
and washed repeatedly with doubly distilled 
water until the wash water reached pH 7. 
The precipitate was filtered, dried at 1 10°C 
for few hours and was analysed 
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complexometrically by EDTA 9 . Densities 
were determined using toluene 10 as a 
solvent and molar volumes were calculated. 

Infrared Absorption Techniques 

Samples used for IR studies were 
washed with acetone and air-dried. Ail the 
bands were recorded on a Perkin-Elmer 
Infrared Spectrophotometer-577 in KBr. A 
few milligrams of the sample was ground 
with two drops of Nujol in an agate mortar. 
About 50 mg of a fine polyethylene powder 
(VESTOLENA 6016 Chem. Werke Huels 
Germany) was added. The resulting paste 
was melted rapidly at about 140°C and 
lightly pressed between two glass plates to 
get slightly wedge shaped film of an 
average thickness of 0. 1 mm. 

X-Ray Diffraction Techniques 

After dying 110°C and then heating 4 
hrs at 950°C, samples obtained that gave 
sharp X-ray diffraction patterns of the pure 
hydroxylapatite phases. The X-ray 
diffraction patters of the samples were 
obtained with Siemens powder 
difractometer with NaCl (TL) counter 
employing Cu-k« (nickel filtered) radiation 
with a scanning speed of l°/min using tube 
voltage of 20 Kv and 24 mA. 

Lattice Constant Measurements 

Cadmium Hydroxylapatite and calcium 
hydroxylapatite and their solid solutions 
with arsenate are hexagonal with two lattice 
constants ao and co. These were determined 
for few suitable samples by measuring the 
diffraction angle, 20 of the three planes, 
(312), (213) and (321). Each sample was 
thoroughly mixed with 25% NaCl 
(recrystallised from HC1) which served as a 
standard so that the observed values of Sin0 
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for the solid solution lines could be directly 
corrected for absorption and instrumental 
errors. The lattice constant of NaCl at 26 °C 
was taken to be 5.6403A 11 . A least squares 
calculation on the corrected values of Sin0 
for the three reflection gave the two 
parameters, ao and Co for each sample. The 
average probable error in unit cell 
parameters is less then +.005 A. 

Results and Discussion 

The results of chemical analyses of the 
sample are given in Table 1. Molecular 
formula for each of the samples was 
assigned on the basis of the results of 
chemical analysis. The molecular formulae 
of the samples were calculated from the 
results of Table 1. Observed values of 
Molar g atomic ratio of Ca/P, Ca/As, Cd/P, 
Cd/As for end members and (Ca + Cd)/ (P 
+ As) were found to be equal to 1.666 
(theoretical 1.67) consistent with the 
formation of homogeneous solid 
solutions 10,11 . This was further supported by 
the closeness of the values of molar 
volumes of the end members and those of 
the intermediate samples which lie within 
the range of end members. 

The formation of solid solutions can be 
confirmed further by the X-ray diffraction 
analysis 12 . Lattice parameters showed unit 
cell contraction / dilation consequent upon 
the introduction of smaller ion Cd 2+ 
(0.97A) and larger ion As 0 4 3 '(2.48A) in 
place of Ca 2+ (0.99A) and P0 4 3 " (2.3 8 A) 
respectively in the apatite lattice. This 
suggests that these are real solid solutions 
and are not mixtures of two components or 
one component with adsorbed material. 

Infrared spectra of mixed hydroxy- 
lapatites were studies in the range 4000-400 
cm' 1 . The results of the IR measurements 


are collected in Table 2. The general 
formula of the apatite considered here is 
Mio(R0 4 ) 5X2 (M=divalent cation like Ca, 
Cd; R0 4 = orthophosphate P0 4 , arsenate 
As 0 4 ; X = monovalent anion, OH). In this 
formula, when M or X is replaced by a 
related ion, the position of the phosphate 
bands in the IR spectra are changed. 

The CaHA and its solid solutions with 
cadmium and arsenate are interesting 
examples of solid state of fairly simple 
molecules of tetrahedral symmetry 13,14 . 
This is the symmetry imposed by the static 
field of the surrounding ions. In addition to 
this field, there is a dynamic interaction 
between the ions, which has also strong 
influence on the internal vibration of P0 4 3 '. 
Under this ideal symmetrical condition only 
two infrared active modes are 
observable 15,16 in regions 1000-1100 cm* 1 
(V3 the P-0 stretch) and 500 - 570cm* 1 (v 4 
and P-0 bending). 

From the spectra of the samples (Table 
2) it is evident that the V3 and v 4 frequencies 
corresponding to P0 4 3 ' are V3 (1075cm" 1 ) 
and v 4 (570 cm" 1 ) for CaHA, V3 (1065 cm" 1 ) 
and v 4 (565 cm" 1 ) for CdHA and 

frequencies corresponding to As0 4 3 ‘ ion are 
V3 (800 cm' 1 ) and v 4 (400 cm" 1 ) and for 
solid solutions containing arsenate, the 
frequencies lie in the above mentioned 
region. The shape of the peaks was also 
effected by the introduction of cadmium 
and arsenate ions into the samples. 

The frequency V3 corresponding to OH 1 " 
ion is observable at 3550 cm" 1 (Cal. 3578 
cm' 1 ) for CaHA. In case of CdHA, the free 
OH 1 " stretching mode appears at 3525 cm" 1 
where as in the case of solid solutions 
containing arsenate, the absorption band of 
OH 1 " ion appeared at lower frequencies. 
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Table 1 

- Chemical analysis of solid solutions fo cadmium - calcium hydroxylapatites containing arsenate Caio- n Cd n 


(PCUVm (As0 4 ) m (OH ) 2 where m = 1, n = 

1,2,3 ■ 

etc. 





SI. 

COMPOUNDS 


Wt.% 


G. 

Molar 

No. 






Atom 

Volume 



Ca 

Cd 

P 

As 

Ratio 


0 

Cag 9 91 (PO 4 ) 5 . 987 (OH) 2 

39.921 

- 

18.487 

- 

1.671 

337.732 

1 

Ca 9 . 004 Cdo .99 1 (PO^.OOS ( As04)q.99 1 (OH ) 2 

32.241 

09.952 

13.850 

6.633 

1.667 

351.984 

2 

Ca7.99iCdi 8 gg(P 04 ) 5.010 (As0 4 )o.989(OH) 2 

26.892 

18.763 

13.029 

6.221 

1.662 

355.516 

3 

Ca6.993Cd 2 .99l(P04)5.00l(As04)o.993(OH ) 2 

22.184 

26.611 

12.260 

5.888 

1.668 

356.900 

4 

Ca5.99 5 Cd4.00l(PO4)4.993(AsO4) 0 .995(OH) 2 

17.978 

33.652 

11.571 

5.577 

1.671 

352.189 

5 

Ca5.ooiCd4.993(P04)5.o03(As0 4 )o.998(OH) 2 

14.220 

39.819 

. 10.994 

5.304 

1.669 

367.062 

6 

Ca 3 .997Cd 6 .006(PO 4 ) 4.991 (As0 4 ) 0 .99o(OH) 2 

10.817 

45.589 

10.439 

5.008 

1.672 

369.301 

7 

Ca 2 .987Cd6.998(P0 4 ) 4 .998 (As0 4 )o.99l(OH) 2 

7.710 

50.662 

09.970 

4.781 

1.671 

370.580 

8 

Ca 2 .00lCd7 97 9 (P0 4 ) 4.995(As0 4 ) 0.987(OH) 2 

4.940 

55.248 

09.550 

4.518 

1.668 

372.131 

9 

Cai.005Cd 8 .975(PO4) 5.005(AsO 4 ) 0.979(OH) 2 

2.377 

59.537 

09.148 

4.328 

1.670 

366.780 

10 

Cd 9 99o(P0 4 ) 5.012(As0 4 )o.986(OH) 2 


63.459 

08.767 

4.172 

1.668 

368.119 

Table 2- 

- X-ray and ir spectra of solid solutions of cadmium - calcium hydroxylapatites containing arsenate. 


Sample 

M f Unit cell volume 

iN0 ‘ 01 Lattice Parameters (A) 

V3a 2 c 

2 


Frequency (cm' 1 ) 


Table-I 







- 

a c c/a 

- 


P0 4 3 ' 

As0 4 3 ' 

OH' 


- 

- 


v 3 

v 4 

v 3 V 4 

v 3 

0 

9.372 6.862 0.732 

521.970 


1075 

572 

— 

3550 

2 

9.272 6.764 0.730 

503.595 


1077 

575 

805 405 

3500 

4 

9.346 6.775 0.725 

512.497 


1075 

575 

808 405 

3475 

6 

9.355 6.755 0.724 

511.969 


1070 

570 

805 402 

3450 

8 

9.327 6.706 0.719 

505.217 


1075 

580 

810 402 

3450 

1° 

9.336 6.677 0.715 

504.004 


1068 

572 

800 400 

3425 
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Abstract 

Chaotic evolution in Froude pendulum driven by 
a periodic forcing has been investigated in detail. 
Control of such chaotic motion has been studied 
by using periodic parametric excitation method. 
(Melnikov’s method has also been employed to 
derive meaningful results. Numerical results have 
been obtained to support the analytical 
calculation. 

(Keywords : nonlinear / chaos / chaos control) 


Introduction 

In the past few years’ extensive studies 
have been done on chaotic dynamics of 
typical nonlinear system. With the results 
of such studies, it has been observed that 
the chaotic motion emerging in the system 
have been suppressed or synchronized .A 
number of techniques are coming out which 
are applicable to control chaos. While the 
dynamics of Duffing equation and Vander 
Pol oscillator are almost exhaustively 
investigated 1 , there have been few studies 
in the system called Froude pendulum. 
Froude pendulum is essentially a pendulum 
mounted on a rotating shaft. The shaft and 
the pendulum are coupled through friction, 
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which is nonlinear. Such system can have 
self-excited oscillation , on account of the 
nature of damping. The objective of the 
present paper is to search a method to 
control chaos in Froude pendulum. We 
have attempted to apply periodic parametric 
excitation procedure here and obtained 
some significant results. 

Chaos control in Froude Pendulum 

The equation governing the dynamics 
of the pendulum can be written as 3 

x + A sin x + cj\x fax 2 - 1) = 0 (1) 

The system represented by equation (1) 
is found to show limit cycle behavior, 
typical of the nonlinear damping term. In 
one study, Thomas and Ambica 4 , used an 
external periodic forcing f costot and by 
writing equation (1) as 

x + A sin x = - q\ x fax 2 -\ ) +f cos «/ (2) 

investigated occurrence of chaotic evolution. 
They have obtained some interesting results 
on this system. To study “smale horseshoes 
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dynamics” of the Froude pendulum, 
Melnikov's technique 5 ' 9 was used. This 
technique is now considered as a powerful 
technique to study nonlinear system. 

In our work we seek to inquire whether 
such chaos occurring during evolution of 
system (2) could be controlled by periodic 
parametric excitation. Thus we investigate 
whether; 

x +a( 1 +cos pi ) sin x=q\ x (<r/ 2 X 2 -l )+/cos cot (3) 

can evolve regularly or not, when A (a 
constant) is replaced by a(l+cos pt) 
(periodic function ) in (2) 

We consider the equation 

x = — <x( 1 +cos pt ) sin x+q\ x (p 2 X 2 -l )+/. cos 
cat 

The system can be equivalently written 
as 

x = v 

v = -a(l+cosp/)sinx + {q\ x(q 2 x 2 ~\) 

+f COS fa)/} (4) 

The undamped and undriven system 
corresponds to 

dx _ ' 

~dt~ V 


— = -a( 1 +cos pt) sin x (5) 

dt 


From (5) we derive the Flamiltonian 
function given by 

v 2 

H(x,y) = — + a(l+cos pt)(\ +cos x) 

For H = 0, we have 

v 2 x 

— = - a(l+cos pt) 2 cos 2 — 

From which, after some calculation, we 
obtain 

x = ± 2 tan' 1 sin h {[(2 kip) sin pt/2) 

v = ± 2 sec h {( 2k / p) sin pt/2) k. cos pt/2 

In this case Melnikov function is given 
by 

oo 

M(t 0 ) = Jv 0 (t-ta)(-q\v ( q 2 v 2 -l ) +/c os 

-oo 

oat)dt 

Evaluating the above integral we get 

M (to) ^ 0 

Numerical Results and Discussion 

We have studied the effect of periodic 
parametric excitation in forced Froude 
Pendulum. 

We have considered the equation 

x +a(l+Xcos pt ) sin x=q \ ± (qp?-\yf cos 
cat 
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Tig. ^-Evolutionary behaviour of Froude pendulum when no periodic parametirc excitation is applied (i.e., Iamda = 0) 
Fig (al. (hi ), (cl) and phase plot represent surface of section, poincare map, time series plot and phase plot respectively. 
The parameter values are a = 1.0, q { = 0.3, q 2 = 0.5, (o = 0.9, /=1.7, p = 1 .8, epsilon = 0.0007. 


<a2) 



1.93106 
1.93106 
1.93105 
1.93105 
1.93104 
1.93104 
1.93103 
1.93103 

HD. 0130. 013^9.01316. 01330. 01316. 013 15 



,'c2 T fflASEFLOT 




Fig. 2-Evolutionary behaviour of Froude pendulum when periodic parametirc excitation is applied (i.e., Iamda = 0.8) 
Fig. (a2, (b2), (c2) and phase plot represent surface of section, poincare map, time series plot and phase plot 
respectively. The parameter values are a = 1.0, q\ - 0.3, q 2 = 0.5, o = 0.9, / =1.7, p - 1.8, epsilon = 0.0007. 
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for numerical calculation. When this effect 
is not into consideration (i.e.. X =0) the 
system evolves chaotically which we can 
observe through surface of section, poincare 
map, time series and phase plots as shown in 
Fig. 1 . Then we start increasing the value of 
X. We observe that for X upto 0.5, system 
evolves chaotically. But for 0.6< X < 1.0, 
chaos in the system is substantially reduced. 
The best result is obtained at X= 0.8. The 
result has been shown in Fig 2. Our 
analytical work also shows .that Melnikov 
integral is zero when X = 1. But when X = 
1.1 the system again evolves chaotically. 
This indicates that chaotic motion in the 
system can be controlled by suitably 
selecting the values of X. 
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Abstract 

Magnetic storms occurring during local midday 
hours are shown to produce much larger decrease 
of the geomagnetic H field at equatorial electrojet 
station, Kodaikanal compared to the decrease of 
H at off-equatorial low latitude station Alibag 
situated in the same longitude sector. These 
storms are associated with the disappearance of 
equatorial Es configuration on the Kodaikanal 
ionograms indicating a reversal of ionospheric 
electric field in the ionosphere during the periods 
of minimum H field. 

The analysis of H data from the chain of stations 
in India shows that the storm time decrease of H is 
almost constant with latitude for local night time 
storms. But for the day time storms the H field 
decreases progressively with decreasing latitude 
from Alibag-Annamalainagar- Trivandrum. This 
clearly indicates the effect of a westward electric 
field imposed over the low latitude ionosphere 
during the main phase of the magnetic storm. 

The ionospheric current contribution identified by 
SqH and the ring current contribution identified 
by Dst -H index are subtracted from the observed 
hourly values of H on the five ID and five IQ days 
of each month of the years 1988-1990 for 
Trivandrum and Alibag. The modified 
disturbance solar daily variation SD (H) showed 
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usual dawn maximum and dusk minimum at both 
the stations. But an additional very distinct 
minimum is observed in AH during the forenoon 
hours at Trivandrum. Difference between SD (H) 
at Trivandrum and Alibag showed largest 
forenoon decrease during equinoctial and least 
effect during June sotstices or during local 
summer months. 

It is concluded, that during the main phase of the 
storm when the ring current is in the developing 
phase, an electric field opposite in direction to the 
atmospheric dynamo electric field is imposed on 
the equatorial ionosphere. 


Introduction 

Egedal 1 has shown an abnormally large 
enhancement of the daily variation of the 
geomagnetic horizontal field, H, over the 
magnetic equator. Chapman 2 suggested this 
phenomenon to be the effect of a band of 
eastward flowing current in the ionosphere 
within ± 3° dip latitudes during the day 
time hours and named it as Equatorial 
Electrojet Current (EEJ). Baker and 
Martyn 3 explained this EEJ as due to 
enhanced electrical conductivities in the 
atmospheric E region within ±3° over the 
magnetic equator. In India equatorial 
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electrojet geomagnetic observatories have 
been operating at Trivandrum (dip 0.6°S) 
since 1957, at Ettayiapuram (dip 2.2°N) 
since 1920, at Kodaikanal (dip 5.0 °N) 
since 1950 at Annamalainagar (dip 7.8 °N) 
since 1957, besides other low latitude 
observatories at Hyderabad (dip 21.8 °N) 
Alibag (dip 25.4°N), Ujjain (dip 33.9°N), 
Jaipur (dip 45.4 °N ) Sabhawala (dip 49.6 
°N) and at Gulmarg (dip 51.6 °N). This 
chain is supplemented by observatories in 
former USSR namely at Tashkent, Ashkha- 
bad, Alma Ata, Karaganda, Navosobirsk 
and others providing a chain of stations 
from - equator to pole within a limited 
longitude sector, unique anywhere in the 
world. Rastogi 4 has described the various 
magnetic and ionospheric effects over a 
wide latitude and altitude regions due to the 
interaction of electric and magnetic fields 
within this relatively narrow region of the 
electrojet. 

Chapman 5 had concluded that the 
magnetic disturbances do not produce any 
abnormal effects on the disturbed, daily 
variation of the equatorial electrojet as the 
quiet day solar daily variations, Sq. It was 
soon observed that the size of the storm 
sudden commencements, SSC, are consi- 
derably enhanced at Huancayo compound 
with those at Cheltanham 6,7 . Rastogi et al. 8 
showed that the latitude variation of SSC in 
H was closely related to the latitude 
variation of the electrojet in any longitude 
sector. Rastogi showed that the SSC at 
Huancayo were associated with simul- 
taneous increase of the electric field in the 
E region of the ionosphere. Rastogi 9 
showed that the Disturbance Daily variation 
of the equatorial electrojet current at 
Trivandrum has a very clear anti-S# shape 
with a minimum around noon, contrary to 
the earlier results of SD variation at 
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Huancayo by Chapman 5 . Later, the study of 
storm time variations of H at the geomag- 
netic observatories in Indo-Russian chain 
showed progressive enhancement at the 
magnetic equator over and above the 
normally expected decrease of Dst ( H) due 
to the ring current effects corrected for 
geomagnetic latitudes. The present paper 
discusses the storm time variations of H at 
Indian chain of stations due to the magnetic 
storms with the main phase occurring 
during the local midnight and midday 
hours. 



Fig. 1- Storm time magneto grams at Kodaikanal and 
Alibag for the storms on 06-15 July 1959 and 
28 November 1959. 
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Fig. 1(a) reproduces the tracings of H 
magnetogram at Kodaikanal (dip 5.0° N) 
and at Alibag (dip 25.4°N) during the SC 
type of storm starting at 1302 hr 75 °E MT 
on 15 July 1959 with its main phase 
minimum H occurring during the local 
midnight hours. The storm has a range of 
785 nT at Kodaikanal but the trace during 
minimum H at Alibag had gone out of the 
paper and so the data was lost. It is seen the 
amplitude of SSC occurring during midday 
hours was 118 nT at KOD and 79 nT at 
Alibag indicating electrojet enhancement of 
day time SSC-H. It is to be noted that the 
short time fluctuations of H are enhanced at 
KOD during the daytime but the night time 
fluctuations in H are of the same magnitude 
at the two stations. The slower component 
of the H variations at the two stations were 
also of the same order. 

Fig. 1(b) shows the H magnetogram 
traces at Kodaikanal and Alibag following 
a SSC at 0452 hr 75°EMT on 28 November 
1959. The amplitude of SSC was 28 nT at 
KOD and 31 nT at Alibag and were of the 
same order as the SSC had occurred during 
the night time when the electrojet was 
absent. The main phase of the storm started 
at 0900 hr and the minimum of the H field 
was observed at 1700 hr LT. The 
fluctuations of H field was much larger at 
Kodaikanal than at Alibag and even the 
general decrease of the H field was much 
larger at Kodaikanal. The storm range of H 
field was 402 nT at Kodaikanal and only 
204 nT at Alibag. The normal midday 
maximums of AH during the noon hours 
observed at Kodaikanal and Alibag on 
magnetically quiet days was replaced by a 
deep minimum at noon hours on 28 
November 1959. A counter electrojet event 
had occurred and the normally observed q- 
type of sporadic E had disappeared during 
the maximum phase of the storm indicating 
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a reversal of the ionospheric electric field 
due to the storm. This indicates that the 
magnetic storms cause the imposition of a 
westward electric field on ionosphere 
during the main phase causing an additional 
decrease of the equatorial magnetic field in 
addition to the decrease due to the 
westward equatorial ring current in the may 
magnetosphere. 



Fig. 2(a)- Storm time variations of the geomagnetic 
horizontal H field at Indian stations 
following a SSC storm at 0109UT on 20 
February 1992. 



Fig. 2(b) Latitudinal variation of the storm time AH = 
Hd - SqH - Dst H at 0030 LT and 1230 LT 
on 21 February 1992. 
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Table 1- 21 February 1992 H storm time 


Station 

Dip 

■ Angle 

0030 

LT 

Dst 168 

Cos im 

1330 

LT 

Dip 

Lat. 

Dst 187 

CosXm 


CosX m 

> 

TRD 

0.6 

-143 

168 

-339 

0.3 

187 

P 

oo 


ETT 

2.2 

-145 

168 

-333 

1.1 

187 

0.6 


KOD 

5.0 

-122 

168 

-267 

2.5 

187 

0.9 


ANN 

7° 8' 

-133 

168 

-229 

3.6 

187 

1.9 

.999 

MUB 

21°49' 

-149 

166 

-181 

11.3 

185 

2.0 

.990 

ABG 

25°23' 

-133 

166 

-168 

13.4 

184 

9.7 

.986 

UJJ 

33.9 

-111 

163 

-142 

18.6 

181 

14.0 

.920 

SAB 

45.9 

-116 

157 

-160 

27.3 

175 

20.9 

.934 


GUL 


Fig. 2(a) shows the storm time 
variations . of the H field at Indian 
geomagnetic observatories following the 
SSC at 0609 LT on 20 February 1992. The 
Dst index indicated significantly positive 
values during midday hours on 20 February 
and the main phase decrease started at 1 500 
hr on 20 Feb. followed by a minimum of 
Dst around the midnight of 20-21 Feb. Dst 
index tried to recover but was followed by 
another energization of the ring current at 
09 hr on 21 Feb. followed by another 
minimum of Dst at 1200 hr (noon) on 21 
Feb. The storm time variation of H at 
Sabhawala (SAB) a station close the Sq 
focus latitude faithfully followed the 
variations of Dst index with minimum of H 
around the midnight and later during the 
noon of 20-21 Feb. 1992. At 
Annamalainagar (ANN), station close to 
the edge of the electrojet belt, the midday 
decrease on 21 Feb. was significantly larger 
than the decrease during midnight of 20-2 1 
Feb. 1992. At Kodaikanal the storm time 
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decrease during the night time was -170 nT 
but during the following noon the storm 
time decrease was -390 nT, even though 
Dst index at these two times were not too 
different. At Trivandrum too the midday 
decrease of the H field was much larger 
than the night time. 

In Fig. 2(b) are shown the latitudinal 
variation of storm time H field at 0000 hr 
and 1200 hr on 21 Feb. 1992 compared 
with the latitudinal variation of the 
corrected Dst index. It is seen that the 
midnight storm time AH showed a small 
decreasing trend with decreasing latitude. 
The midday storm time AH showed a small 
decrease from Ujjain (UJJ) to Alibag 
(ABG) and Hyderabad (HYB) all outside 
the electrojet belt. There was a strong 
decreasing trend of storm time AH 
thereafter from Hyderabad - Annamalai- 
nagar-Kodaikanal - Ettayiapuram to Tri- 
vandrum, close to the equator. The 
observed decrease at TRD was -340 nT 
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compared to the Dst index of only -170 nT. 
This clearly indicates that a strong 
westward electric field was imposed on the 
eqauatorial electrojet region during the 
main phase of the storm. 



Fig. 3(a) Modified Disturbance Daily variations of 
the horizontal magnetic field H, (SDH)* at 
equatorial electrojet station, Kodaikanal and 
at off equator low latitude station, Alibag 
average for there seasons of the years 1988 
to 1990. 


1988-1990 
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Fig. 3(b) The difference between (SDH) * variations 
at Kodaikanal minus the same at Alibag for 
different seasons of the years 1988-1990. 

Next, from the monthly data tables of 
the observed hourly mean H field at 
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Trivandrum and Alibag, another table was 
generated giving the hourly values of the 
observed H field minus the mean H field on 
five international quiet days (i.e. SqH) 
minus the corresponding hourly mean value 
of Dst index. This modified (AH)* indicates 
the remenant H field at three stations 
devoid of the effects due to normal quiet 
day electrojet current and due to the 
disturbed day ring current. Next, mean 
daily variation of this modified (AH)* on 
five 1 Q and five 1 D days were computed. 
Next, mean of the variation on ID minus 
the same on IQ days were computed giving 
the modified Disturbance Daily variation 
(SOD)* of the H field at these stations. Fig. 
3(a) shows the modified (SDH)* averaged 
for the three seasons and for the whole year 
for the period January 1988 to December 
1990. It is to be noted that the yearly 
average (SDH)* at ABG and TRD show 
comparable values during the hours midday 
to next sunrise, with the usual dawn 
maximum and dusk minimum. During the 
forenoon period from 0600-1200 LT 
(SDH)* at ABG shows slow decrease from 
dawn maximum to dusk minimum. But 
(SDH)* at TRD shows a very distinct 
minimum during this period of the day. 
Examining the seasonal effect on this, it is 
seen that the forenoon decrease of (SDH)* 
at TRD is most pronounced during 
equinoctial months, less during D months 
(local winter) and least during J months. 


In Fig. 3(b) are shown the difference of 
(SOH)* variations at TRD minus the same 
at ABG, thus removing the common cause 
of disturbance daily variation at the time 
stations namely the dawn-duck magneto- 
spheric electric field effect. A very clear 
forenoon decrease of the difference curve is 
seen during the equinoctial and December 
months. 
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It is suggested that this additional effect 
is due to the immediate transfer of the polar 
electric field during the phase of the 
magnetic storms when the equatorial ring 
current is being emergized. 
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ACADEMY NEWS 


(I) AWARDS 

1. Prof. R.S. Sirohi, Ph.D., f.n,a.e„f.n.a.Sc., Director, Indian Institute of Technology, Delhi 
has been selected for the Hari Om Ashram Trust Award entitled "Sir .C.V. Raman 
Award: Physical Sciences", for the year 2002 by the University Grants Commission. 
This award includes a sum of Rs. 50,000/= and a citation. 

2. Dr. V.P. Sharma, d. P hi/, d.Sc.,f.n.a.,f.a.Sc.,f.n.a.Sc.,f.a.m.s.,f.ra.s, formerly Additional 
Director-General, ICMR & Director, Malaria Research Centre, Delhi has been conferred 
24 th G.P. Chatterjee Memorial Award (2003-2004) by Indian Science Congress 
Association, Kolkata. 

3. Dr. Jayant Vishnu Narlikar, Ph. d. (C antab.), Sc. d. (C antab.), f.n.a.,f.a.Sc.,f.n.a.Sc., f.t.w.a.s., 
Director, Inter-University Centre for Astronomy and Astrophysics, Pune, has been 
conferred Padma Vibhushan by the President of India. 

(II) PROGRAMME OF SCIENCE COMMUNICATION ACTIVITIES OF THE ACADEMY 


Popular Lecture January 27, 2004 by Prof. Govindji at Allahabad. 

Science Extension Lectures February 10-20, 2004, at Allahabad 


Outstation Science Extension 
Lectures 


Local (Allahabad) Level Science 
Communication Activities 


February 10-15, 2004 at Satna & 
adjoining areas of M.P. 

February 15-20, 2004 at Jaunpur, U.P. 

Science Fiction and creative contest on 
17 th February, 2004 


State Level Science Contests 
(to be held at Allahabad) 

State Level Essay Contest 
State Level Teachers Workshop 
National Science Day 


Science Quiz February 24, 2004 
Science Debate February 25, 2004 
Science Oration February 26, 2004 
Science Exhibition February 27, 2004 

For Degree Colleges of U.P. on February 27, 2004 

From February 25-27, 2004 at Allahabad. 

February 28, 2004 being jointly organized with 
Harish-Chandra Research Institute, Allahabad. 
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FORTHCOMING SYMPOSIA/SEMINARS/MISCELLANEOUS ANNOUNCEMENTS 


1. G.D. BIRLA AWARD FOR SCIENTIFIC RESEARCH 

Last date for sending Nomination 31 st March 2004 

For details please contact 

K.K. BIRLA FOUNDATION 
Hindustan Times House, 

10 th Floor, 18-20, Kasturba Gandhi Marg, 

New Delhi - 1 10 001 

Phone : (011) 23718282, 23317735 


2 - THE DISABLED CHILDREN ASSOCIATION AWARD FOR SCIENTIFIC RESEARCH 
(20041 

Last date for sending Nomination 30 ,h November, 2004 

For details please contact - 

The Award Secretariat 
Disabled Children Association, 

P.O. Box -8557, 

Zip Code : Riyadh 11492 
Saudi Arabia. 
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